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XII. On Stokes's Current Function, 

By E. A. Sampson, B.A,, Fellow of St. John^s College, Cambridge, 

Communicated by Professor Geeenhill, F,R,S, 

Received ."N"ovember 24, — Read December 11, 1890. 

In Maxwell's ' Electricity and Magnetism ' ^ a view is put forward, in accordance 
with which we may regard any irrotational motion in a perfect liquid, for which the 
velocity potential is a solid zonal harmonic, as due to the juxtaposition at the origin, 
and upon the axis of symmetry, of sinks and sources. 

But, in a liquid, any irrotational motion which is symmetrical with respect to an 
axis gives a velocity potential which may be expressed as a sum of a series of solid 
zonal harmonics, their common axis being the axis of symmetry, and their origin 
arbitrary, provided it is excluded from the region to which the expressions apply. 
The position of the origin upon the axis is arbitrary, since, by a transference formula, 
we may pass from one origin to another. 

Let us now consider the system formed by a line source and a line sink, of equal 
strengths, extending along the axis from an arbitrary origin to infinity in opposite 
directions. Such a system I shall call an extended doublet, of strength m, where m is 
the strength per unit length of that part which lies on the positive side of the origin. 




By the superposition of two extended doublets of equal but opposite strengths we 
can produce a sink or a source upon the axis. Hence, in a liquid, any irrotational 
motion which is symmetrical with respect to an axis may be produced by superposition 
of extended doublets, whose origins depart but little from an arbitrary point on the 
axis of symmetry. 

Now, for an extended doublet of strength m, Stokes's Current Function xjj, for 
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any point distant r from the origin, is — 2mn For let ^ be the distance of the origin 
of the doublet from the origin of coordinates, and let xfj (m, Q be the value of Stokes's 
Current Function for any point (m^ z). 




Then if S\p be the Current Function for a source of strength 2m BC at the point C of 

the axis we get 

1 fJ 



1 J^ 
us rdd 



Sxjj 



2mS^ 



Therefore 



cie 



Bifi -= 2m — 8{, 



= 2w sin ^8^5 

whence 

Sxjj = — 2m8^ cos 0, 
the constant of integration being zero. 

But 

Sxjj = xjj {m, ^) + i/; {— m, I + S^), 

and clearly, 

i/f ( — m, C + SQ == — i|f {ill, C + SC). 
Hence 



tfj (m, 



^i^hO + j^H 



= — 2m8C cos 6, 



2mBC 



z-K 



^W' + (^ - m 



Therefore 



^-? 






91 ? 



and 



V' 



i:mr 
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where 

disregarding a constant. 

Thus, if 

m=f{C)dC, 

we may, by properly choosing the function /, write 

rP= f(0^i^^-^{z-Cf]dZ (2), 



00 



where \p is the current function for any irrotational motion in a Hquid, symmetrical 
about the axis of z. 
Again, if 



ineretore 
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dr z -- ^ 
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- r' 


dz r ' 


dh' 


1 m^ 


^V 1 (z tf 


dm^ ^' 




dz^ r r^ 


dh^ 
dm'^ 


d^r 2 
dz^ r 
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nr^ + (z ^f 
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> 




1 


dr 




m 


dm^ 



and the expression (1), and consequently also (2) satisfies the differential equation 

dm'^ dz^ m dm ^^^^ 

or, as I shall write it, Dt/r = 0. 

When the motion is rotational, (3) no longer holds. In fact, as is well known, we 

have under all circumstances 

1 

— J)\h = — 2ft) , 

w ' 

where &> is the resultant spin at the point {m, z). 

Thus, if there is spin in the fluid, (3) is replaced by 

am^ . dz"^ m dm ^ ^ 

3 M 2 
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Again^ if v^ stand for the operator 



Cv . tv X ih . JL 

dm'^' dz^ m dm m 



^ d(j)^ 



(f) being the azimuthal angle about the axis of symmetry ^ it may be seen at once that 



B^ 



Consequently, (3a) may be written 



sm <p tu ^ 



B • d * 4 « 9 



(4). 



3 



'^siiKi 



tjy 



2o> sin <ji 



d o 



s a 



. (3a). 



Consequently^ 



^ = -^0 + 



m 



27r sin ^ 



' sin (^^ft)' (^^' dy^ dz' 



tJ tj el 



r 



where ijjo is a solution of (3), or, xjj consists of a solution of (3), together with 

— ^r~ X the potential at the point considered of a distribution of mass of density 
27rsm(^ ^ 

at any point, sin<^ X the spin at that point. This result is given by Basset, 
^Hydrodynamics,' vol. 2, § 306. 

I give one other general result. Since 



1 



CO 



2 



— Dxjj 

'US 



* a u a 



s u 



. . (3a), 



the circulation in any evanescible circuit drawn in a meridional plane is 



J », 



" Dxjj dm dz . 



P?'«9eo0a3 



(5) 



where the integration extends over the area embraced by the circuit. 

This result enables us to transform Di// readily from cylindrical to other systems of 
coordinates. For instance, consider polar coordinates, r, 9, and let us find the circula- 
tion in a small rectangle bounded by r^ r + dr, 9, 6 "\- d9. 

Let the velocities in the direction of r, and perpendicular to it, be R, @. 




ME. R. A. SAMPSON ON STOKES'S CURRENT FUNCTION. 



453 



Then the circulation in this circuit is 



Il(ir + 



dr 



®r + ~- (@r) dr 

'd% 



de 



K+^f'' 



dr — %T dd 



r dr dd 



dr r rdd 



Now 



® = — 






1 d'\^ 



T sin d dr 

1 d'\^ 
r^ sin 6 dd 



Thus the expression in square brackets is 



or, 



T sin 6 



Biji 



d^ '^ r^\dd^ "" c^6^ 



^^'^ sin ^ ^ / 1 d^fr' 

'" r 



dr^ 

d^'xjr 



immtmm 



r^ dd \ sin 6 dO 

1 — fA? d^-s^ 



■ (6), 



if \h stands for cos 6, 

Other applications will be found later. 

Reverting now to the expression (1), it will be seen that the direct distance of any 
point from a point on the axis of symmetry, plays the same part in the theory of 
Stokes's Current Function that is played by its reciprocal in the theory of the potential 
function belonging to symmetrical distributions of matter. 

Thus, if <>% 0, r, Q be the coordinates of a point upon the axis, and of any other 
point, the distance between these points, s/{r^ — 2r^Tcos d + r^), may be developed 
in a converged series, say 



M, = 00 

t 



m 



n ~ CO 



fji n 



r, 



n- 



- I^^ (cos &) or 



I 



S ~ ;;;>-3T I» (cos ^) 

71 = ' 



according as r^ is greater or less than r, \ (cos 6) being a certain function of 6, and 
we see from (6) that 



{l-l-')'^ + n{n-l)U,.)^0 . 



(7) 



Now it is evident from the analogue of Zonal Harmonics, that it is proper to discuss 
the function L (cos 6), and other solutions of (7) before considering the applications of 
Stokes's Current Function to the motion of liquids. It is with this discussion that 
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the first three chapters are occupied^ and, as might be expected, the theory closely 
resembles that of Spherical Harmonics. I have accordingly made free use of the order 
and methods adopted by Heine, in his ^ Handbuch d. KugelfuDctionen/ more especially 
in Chapters I. and II., where the necessary changes were slight. Moreover, the 
functions I deal with have themselves been discussed by Heine, on a different method, 
and most of the expressions which I find in the following pages are given by him. 
Full references to these are given on p. 461. 

The idea of developing the solutions of D\// === in a manner more or less analogous 
to that employed with regard to Laplace's equation, appears to have been first used 
by 0. E. Meyeb,,"^^ who obtains the equation (7), shows that the functions contain 1 — /x^ 
as a factor, and that they obey (22a), Chapter IL An expression which shows the 
relation of the functions to Zonal Harmonics was given by Mr. BuTCHER,t and 
functions of fractional order have been used by Mr, HiCKS.J in connection with his 
researches on the theory of the motion of vortex rings. The fuller account of 
such functions, which is found in the following pages, may be of interest in relation to 
these ; for example, I would refer to p. 44. 

The applications to Hydrodynamics, which I here give, are of mathematical interest 
rather than physical. They are chiefly in connection with the motion of viscous 
liquids. In ' Crelle-Borchardt,' vol. 81, Oberbeck has given the velocities produced 
in an infinite viscous liquid by the steady motion of an ellipsoid through it, in the 
direction of one of its axes, and from these Mr. Herman § has found the equation 
of a family of surfaces containing the stream-lines relative to the ellipsoid. In 
Chapter VI., Stokes's current function is obtained by a direct process for the flux 
of a viscous liquid past a spheroid, and it is shown that the result differs only by a 
constant multiple from the particular case of Mr. Herman's integral. 

Some minor applications are also given ; namely, the solutions are obtained for flux 
past an approximate sphere, and past an approximate spheroid. The solution is also 
obtained for flux through a hyperboloid of one sheet, where it appears that the stream 
surfaces are hyperboloids of the confocal system. A particular case is that of flux 
through a circular hole in a wall, and this is interesting, because we see that by 
supposing internal friction to take place in the liquid, we find an expression which 
gives zero velocity at the sharp edge, and thus avoids the difficulty which is always 
present in the solution of such problems, on the supposition that the liquid is perfect. 
A comparison may be instituted between this problem and that of the effect of a 
disturbing periodic force upon a dynamical system capable of vibrating alone with a 
period equal to that of the force. It is well known that the amplitude of the 
vibration induced appears infinite, if we totally disregard friction, and this difficulty 

* ' Ceelle,' vol. 73. 

t 'London Matli. Soc. Proc.,' vol. 8, see p. 49. 

X *Pliil. Trans.,' 1884; 1885. 

§ ' Quart. Journ. Math.,' 1889 (No. 92). 
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is met by the fact that the (lamping effect of even slight friction is increased consider- 
ably by high velocities. Now a viscous liquid can move irrotationally, and if there 
were no friction at the boundaries this is the class of motion it would take in cases 
of flux past or through obstacles. But if the obstacle terminated in a sharp edge, 
this would make the velocity there infinite, and the friction, however inconsiderable 
elsewhere, would here become of account. The boundary conditions which were 
necessary for the existence of irrotational motion throughout the liquid would no 
longer apply, and the whole character of the solution would be changed. 

This would, at any rate, seem to apply to cases in which the whole motion is slow, 
and when, consequently, the boundary conditions which must hold are pretty well 
understood. 

The paper concludes with an attempt to discuss the flux past a spheroid or through 
a hyperboloid at whose boundary there may be slipping. The current function is not 
obtained, all that appears being that it probably differs from the parallel case of the 
sphere in being far more complicated than when there is no slipping. From this we 
except the case of flux through a circular hole in a plane wall, when the solution for 
no slipping satisfies the new conditions. 



Chapter I.— Various Forms of ln{oo). 

Where the method of development in the following pages does not necessarily differ 
from that of Heine, I confine myself to a statement of result, with an indication of the 
method, and a reference to the corresponding section in his ' Handbuch,' vol. 1. In 
one place the reference is to Ferrers'^ ' Spherical Harmonics.' 

Let R denote ^(1 - 2aaj + x% 

Then writing 

we find 



•I;?. V / — 



1 . 3 ... 2^ — 3 
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the series ending with x^ or od 
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n even] 
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2' 2 
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2 ' ^'^ "^ 2 
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n odd] 



^ 



-». 
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'n 
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2/ 
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-^ jO_ /y> 

2? *2? •^^ 



^ 



• (9^) 



y • ■ (9b) 



There are two exceptional cases, viz., 

The first eleven functions are given in Chapter IV., Table II. [Cf. Heine, § 4.] 
If 






and therefore 



CCS 0^ and ^ ^^ x -- \^{x^ —. 1) == ef^ 
2 cos nl9 = t + i'^'\ 



such a sign being given to the square root as shall make mod. ^ :|> 1, and therefore 
such a sign to id as shall make its real part, if not zero, negative, we have 



T /^\ -^ ('^ I) 



'^ 



2lJ(n)U{-^i) 



rj-^i-h 



n. 



n + 1, e) 



n i^l ^. 3.\ 



!>■ 



(10) 



n (?2 -— I) 

iTo^)n(--i) 



2 cos nd — 



2^ 



2 (2^ -- 3) 



« 2 cos (^^ -' 2) 



2^^ (2^ - 2) . 1 



2 cos (^1 — 4) ^ 



^ ' r-— rr . 2 COS 0^ - 6) ^ . , . 



2 . 4, 6 (2?^ -^ 3) (2^i -- 5) (2?^ - 1) 



(10a) 



n being in each case a positive integer, [HeinEj § 6/ 
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e shall show (p. 471) 

T (t\ 



■n /•! 



TT J^i 



^t''^^l - uf (u -- i'^f du ..... (11). 



This form is equivalent to (10) when t^ is a positive integer. 
Putting 



ay =1 -— ^/(l •— 2iax 4" ^^) 5 



so that 



a 



y= a? + -(/-- 1); 



whence^ by the help of Lagrahgb's theorem, 



I4-) = ^ (ly " (S"-r • [Hbinb, § 7] . (12). 



The distance between two points whose coordinates are m^ z; 0^ Cl is 



where 






y^ {tST^ + (^ ■"" Cf} 

zjr , h = t,jr , r = ^/(ot^ + z^) ; 



whence^ by Taylor^s theorem, 



mfl — 1 cf**^" 

I„ («) = (- l)»-i -^ . — [Fberehs, § 5] 



n dz 



• (1^)- 



Heine shews, §§ 8, 9, that 



J ft 



(i'j; 



TT 



qocX — 1 — a COS V\/{^^ — 1) 



y/{l — 2aa? + 54^) 



provided mod. a\/x^ — 1 < mod. (ax — 1), a result which is secured by taking a 
small enousfh. 
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1 -— «£» — s/i^ — - 2aaJ + a^) 
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tlT 
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TT 
1 ''^ 



V 



1 — a^ -j- 



1 — 2£«-:» -4- a^ 



«!«? 



1 -- a\/af' — 1 . 



COS'^ 



TT^ 







diy . 



^3 (^l ^ ^$ ^ Xs/x^ — 1 . COS 7]) — a COS ^-v/^^ "~ 1 

a (a? — COS 7j\/m^ — 1) — 1 



Let each expression be developed in ascending powers of a^ and equate coefficients 

of OL^K 

MDCCCXCI.— A. 3 .N" 
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Then 

I^{x)==^-^ (x^ --- l)8in^ 7] {x -- COB ri\/x^ ■— lY''^dy] , , . {14), 



where 



-J {x^ — 1) sin^ (f){x + cos ^\/^^ — 1)'"^"^^ . . (14a), 



1 =z (cc -^ cos 7;\/x^ — 1) (x + COS ipx/x^ ~~ 1) [Heine, § 9]. 



The form (14a) is inadmissible when (r is a pure imaginary ; for suppose the path of 
integration of rj to include real values only ; then ^ is complex and discontinuous, 
passing at a step from 0+iooto7r+^^- 

The forms (14) and (14a) may be derived from (11) by substituting 

It/^ z=z X— COS 7) \/x^ — 1« 

From (14) we find a form analogous to one of Mehler's forms for the Zonal 

Harmonic. 

Write 

X — cos J) \/{x^ — 1) = e\ 

where 6 will in general be complex, and we shall suppose its imaginary part to lie 

between i ^^r. 

Then 

sin 7) ^{x^ -- 1) = y(-- 1 + 2xe' - e'') 

if such a sign is given to the square root at the right as shall make its real part of 

like sign into ^{x^ — 1). 

Thus 

1 



I^{x) = -j 'd9 y/{2 {x - cosh e)] e(^^-^)' .„...« (15), 



0n 



where 



X ^ \/{x^ -»- 1) — 6^" 
X + s/{x^ — 1) = 6^\ 



Thus if ;r = cosh r9^, and the imaginary part of ^ lie between i ^7^, 

and 

I^^ {x) z::^ " dd ^^{2 (x — cosh 0)} [ cosh (n --- i) ^ + sinh {n — ^) 6] 

s= — " dd ^^{2 [x -^ cosh 6)] cosh (n — J) 6^ . . . , (15a), 
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where sinh (cosh~^ .'r) = ■— y^ (o;^ — 1), a,nd the imagiTiary part of cosh""^ a? lies 
between i ^^r. 
Write ^^ for $ 

I.^^ (.x) ~ ™-« c/(j^ \/{2 (cos (f} —- .q:;)} cos (n — i) (^ , . . (ISb), 

where sin (cos""^ a?) = ^{1 — a^^), and the real part of cos""^ x lies between ± tt. 

Another form may be obtained when x is real and lies between ^ 1, and n is a 
positive integer, provided we make the same assumption that underlies the ordinary 
process of finding Dirichlet's forms for the Zonal Harmonic, viz. :— 

I^^ (x) = ~ sin {n — -1) ^ ^{2 (x — cos (j))} clef). [Heine, §11] . (15c). 

''^J COS~' X 



We have (p. 453), 



{l-x^)'^ + n{n-l)I4x)=:0. 



From this equation we obtain the forms (9), (9a), (9b), (10), directly. 
1st. Put 

1 

-■ « 

then 






if 



f. 1 • 
this 



{m — n) {m + ^ — 1) := ; 



is a solution. 
2nd. Put 






x^ = t 



'(l-«)S + f^-l')f + -<^'l. = (7bV 



of which 



is a solution. 
3rd. Put 



dt^ ' \2 2 dt 



"1—^ / ih ih JL JL 

■^' ~2' 2~2' 2' ^ 



3 N 2 
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of which 



^> ~2 + 2' 2' 2' ^ 



is a solution. 
4th. Put 

Then 



if 



<(l-«)|f+[m + f-(m + i)^]| + ™.= .... (7d) 



(m + n) (m — ^ + 1) = 0, 



Of this 



is a solution. 

From the differential equation (7), we may form a conception of I^«(^), where n is 
not restricted to integral values ; the discussion is found in Chapter III., where 
reasons are given for choosing (H), or its derived and equally general forms (14)^ 
(15)^ (15a)3 or (15b) as the definition of ln{x). These are identical with (9), (10), 
when ^ is a positive integer, and not otherwise. 

I^(dr l) = 0. This is obvious from (11), When n is a positive integer all the 
roots of In{oc) = are real, different, and lie between ±: 1, and excepting zero, occur 
in pairs of the form ^a. This may be seen from (12). [Heine, § 7. 

The sequence equation 

(n + 1) I^^iix) - {2n - 1) xl,, (oc) + (n - 2) l,^^{x) =0 . . (16) 

connects three consecutive functions ; this follows readily from (15b) ; n is not neces- 
sarily an integer. 

We also find — a result required hereafter— 

X^ i;, {x) = S,, L + 3 (x) + en lu {X) + L I..--2 (^) . . . . (l 7), 

where 

^ (n -f 1) (w + 2) 



(2n + 1) {2n — 3) 
(n — 2) (n — 3) ^ 

exceptions, 

We have for arbitrary n 
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Whence, comparing with (15b), 






• • 



(18), 



Other relations, of use in transformations, may be introduced here 



{2n— 1)1^ = P,, — P^_.2 



• « 



~n 



-P, + 2a;P«- 



n^2 



• • • • * IX(7la 



(20), 



• » • • » 



dl 



(1 — x^) ~r-' z= {n—l)x l„ — {n + 1) I„ + , 



dx 



= — nx I„ + {n — 2) I„_j 



« * • c 



(21). 



These hold without restriction upon n. 
We see from (18) that 



3 



1 .2 . , . n 



^^i^'^-"^)(a^)^ 



where ^^^'^"^^ {x) is defined and discussed by Heine (Chapter III., §§30-33, and 
Chapter IV.), who obtains equivalents of the following results :— (9), (10), in § 51 ; 
(12) in §31 ; (14), (14a) in §50. Also the following, which are proved below:— 
(28), (32) in § 51 ; (33) and equivalents of (34), (35) in § 52 ; (22), (23) in § 62. 

These might, therefore, have been assumed and the theory based upon them, but it 
seemed preferable to indicate briefly the direct development. 



Chapter II.— Developments in Terms of I {x). 



From equation (7) we find 

(m — 7^) (m + ^ — 1) 
Whence 



^m \p^) ^n y^j 7 



X' 






dx 



+ 1 



(22). 



+ 



Iw ('-'^7 1% (^) 7 ^ 



-1 



X' 



. . (22a), 



when m ■=[=■ n^ and neither of them is or 1. 

The limiting value of this integral when m = r?, is 



f 









dV 

an 



+ 1 



-i 



2?^ 



dn 



P^«x 



+ 1 



-1 



Using (15b), 



dl 



n 



dn 



2 



IT 



C cos ■*" 1 a? 







(j> . sin {n — i) ^ 'x/ (2 (cos (^ — - oc)] c/^. 
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This vanishes at the upper Umit, and at the lower gives us 

~ \ (f> ^ sin {n — -|-)^ , 2 cos - . d<j). 
Integrathig, we get the required expression 



2n - 1 
and when n is a positive integer 



2 cos ?i7r 2 sin 7^7^ . (27-^ •— 1) 



_n (n — 1) 



n^ (71 — 1)2 






a • 



« • « 



(23). 



Every member of the set of functions, \ ix), \ {pc), . . . I,, (cc) . . . . vanishes with 
X '\- \ and X — \ I besides, \n {x) has n — 2 real factors lying between these two ; 
and any two members of the set obey (22a). Now, if we have another function, ^ {x)^ 
also vanishing with x '\- \ and. x — 1, it has been shown*"'' that we can find a linear 

00 

function of the Ts, SAJ^^((:c), which shall be equal to <^ (^) for all values of x 

2 

between + 1 and — 1, and 



A 



n 



+ 1 / r + 1 

dx , (f){x) . I,, (^)/(l — x^) I \ dx , I,, (x) \ {x)l(l — ^^ 
-1 



'/)-> 



x^). 



Or, 



n — 00 



1/ \ ^ n {n - I) {2n — I) f + i , ^ I„ {y) I„ (x) 



e a 



'/^ = 2 



r 



. (24) 



where 



(^(±1)==0. 



an 



d 



i ^ «!^ ^ 



It is easy to prove that this development is unique ; I shall now show when it is 
convergent. 

From the sequence equation, (16), we find 

{n+ l)n{n-' l)\l,,^^{x)l,,{y) '-ln{x)ln+i{y)\ 

— n{n'- 1) {n — 2) [I,, {x) \^^ (y) - I,^^ (x) I,, {y)l 
= {x — y).n{n'- 1) {2n — 1) l^ (x) 1^ (y) 

from n= 2, onwards ; and, therefore, the sum of the series (24) ton terms is 



dy . <i> (y) . (n + l)n{n-l) 

-1 



(x ^ ^) (1 ~ y^) 



**Liouville/ vol. 1. The method is quoted, substantially, by Lord RAYLEian, * Theory of Sound,' 
vol. 1, eh. 6. 
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Similarly, it may be shown that the sum of n terms of the development of f{or) in 
Zonal Harmonics is 

AjiA. it is shown by Heine (vol 1, 1 1 19), that the limit of this, when n is infinite is 
the mean ^ {f{x + 0) +/(^ -— 0)}, provided f{x) is not a function with an infinite 
number of maxima and minima. 

Moreover, (Heine, vol. 1, § 40), when n is very great 



x^i \Xj 






and it will be shown, (p. 476), when n is very great 

I (^) _ rid- ri' 

where ^ = x — ^/{x^ — 1), as on p. 456. 

Therefore, the two integrals above tend respectively to the same limits, as 



1 



1 ri/.A^^'')"'^^-^')*./^ 



'%) . ^-S^ ! ^S ^. 



where yj =i y --- ^/{y^ — 1). 

The limit of the second we know, and if we write 

the second reduces to the first. Therefore the limit of the first is the mean value of 
^ {x) at the point x, provided ^ {x)j^{l — x^) is a finite function of x, which does not 
possess an infinite number of maxima and minima. 
We find the function 

^ "^ 2 "^ ,?i ^' ^^^ ~~2~ ' ' (t^Tl^"-l)7^T 2r^" 3^ . . . (^ + 1) 

I ' v" T / \ i^l±i (^ •" ^) 0^ - 3) . . . (^ - 2r + 1) . 

+ ^^^ l2,+i (^; 2 • ^^ _^ ,^^,^ (^ + 2r - 2) . . . 0^ + 2) ' ^ ^^''^ 
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is one whicli is equal to q6^^ when x is positive and to zero when x is negative, for all 
values of n. 

When n is a positive integer 



-L. • ^ • S • iV I /■ y-^ __, 



X a O . , , ^'??/ "~~ X 



n 



x) + (2^^ -"-^ 5 



ZiYh '^'^ X 



«-3 



(a;) 



L- 



, ,^ „, (%% - 1) (2% - 3) TT 



(25a), 



whether x is positive or negative ; this eqxiation is in fact an identity. 
And if 

F {x) = Cq + c{x + CgOs*-^ + . = 
we may write 

F {x) = 6oIo (^) + ^ili (^) + %J-3 (^) + » . M 
where 



& 



n 






. 0^ + 1) (% + 2) 



+ 



2 (27^ 4- 1) 
(^ + 1) (^^ + 2) (?^ + 3) (^s% + 4) 

T747(2^r+^ 



'/^ + 4 ' 



. , . IZoBI 



from fi = 2, onwards, while 



We have 



6o + 6i =: F (1) ; - 6, 



cos mC7 = cos -;.-' COS Wl - 



^1 



( — 1), [Heine, § 16. J 



u ] + sm -77- sm m[-r — u 



sm mu = sm — cos m ( ~ — u 



7)l7r a / TT 

cos -;;" sm M . ^ 

2 



and if cos 6 = x, we find 



TT 



cos m 



e 



mir 



COS 



2 



Xq \X} *~~" o 



m' 



m' 



V 



I2 (^) + ^ 



m^ . m^ 



03 



7?i 



3 







33 



I4 ('^) 



TT 



sm m 



^ 



sm 



mTT 



i^ 



\ {x) - 5 



m 



3 



X3 



3 ^ 93 



m' 



I3 (a:;) + 9 



m^ — .1''^ . 7)1? — 3' 



m 



2 _ 2^ m^ ^- 42 



x) 



. (26). 



When m is a positive integer 



COS m$ = 
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2 :. 4 . . . . , 2m — 2 



J..O • • • £J i/li 



(2m — 1) I», (cos 0) 



+ 2:T:2^-¥:i^r3i(2^'^-9)i»-4(coB^)+...J . (26a). 

Where m is even, m^ occurs in the numerator ; where m is odd^ m is absent from 
the denominator. [HeinEj § 19. 
We have, by (19) 

V,,{x) = {2n — 1) ln{oo) + P^^-s (^) = (2ti — 1) L (.^) + (2n — 5) I,,„g (cc) + » . . 
Hence, if we can effect the development 

F (x) = CqPq (x) + e{Pj (x) + <^2^2 (^) + • • • ,^ 
we can write 

F (x) = &oIo(^) + ^Ji (^) + Ms (^^) + . » . , 
where 

h,, = {2n — 1) (cn + C;,4.3 +...). ... . . . (27). 

If we develop (I — y^)l{y "^ ^) hi ascending powers of x^ and then transform by 
(25b), we get 

J 1 1 1 1 Qf/y* ?t> =: 00 

— --—^ zz: — a; + S ^^(tl — 1) (2n — 1) H^, (y) I,^ (i^), 

where 

TT / \ n (^ — 2)11 (—1) ■ ^ . ., ^^ fn "-In . , ^\ ,^^. 

H42/) = - ^Vll^- V^" ^ r¥ ' 2' ''+*' ^") • ^28). 

The convergency of this development is discussed on p. 476. 

In future references to the expression H,, (y) we shall not in general suppose n to 
be a positive integer. On comparing it with the expression (9) for I,, (a?), it is seen 
that the former may be derived from the latter by changing x into y^ n into — n + 1, 
and multiplying by a constant factor. 

Since a similar change transforms 

(1 -oc^)'^, + n{n^ l)z = into {1 - y^)~ + n{n - l)z - 0, 

we see that H,, (x) is another solution of (7), a result at which we have already 
arrived on p. 459, 

There are two independent solutions to the equation (7), whatever the value of the 

MDCCCXCI. — A. 3 O 
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independent variable may be, but (28) ceases to be a permissible expression for one 
of them when mod. ^/^ < 1. 

For mod. 2/^=1, the series is convergentj since (ti — l)/2 + %^2 ■— (^ + |) is 
negative, and we have, by the well-known expression for the sum of a hyper- 
geometric series of which the fourth element is unity, 



^n \^) 



n (n --" 1) 



e f» a • « 



8 i 



. (29). 



The function H^ (m) is of odd degree in x when n is even, and of even degree when 
n is odd ; it vanishes when x is infinite. For n = or 1 the given expression fails, 
but if H is still a solution of (7), different from the corresponding I, we may put 



XI A \Xj "— X J JUi \Xj — -"■ X ■ « 



♦ « • » • 



(30). 



Certain functions which are not finite for all values of x between ± 1, may be 
developed in a series involving H (x) in place of I (x). 



^=t n(n ^ 1) (2n « 1) H, (x) -^^^^ = 2 i2n - 1) H, (x) ^^^ • 



2 3 



differentiate n — 2 times with respect to y, and then put 3/ == 0, and we find 






n + l 



A. ,^ , , t (1/ """ ^ 



1 • rt 



• . • 



we may re-wi*jte it 



where 



'n 



2-^—1 
{2n — l)lln{x) — (2n + 3) ^^^-"y" H,, ^, ^ («) 

_1_ /Q/T) JL ^\ ^'^ — X.2fl + 1 TT / X 



« « # 



= 62H3 (O^) + Sglig (^) + * « • ? 



X • O » « » -ij^ "~~" o 

1.2 . . .^1 — 2 



f?- — 2.f^/ «— 3 



Gfi^ 



n^l 



/j . ^'?1!' 



iy ^}i B»* 3 



o 






n 



2.% — 3.^ — 4,n — 5 



2.4. 2^^ --3.2^1 — 5 



/% 



#-5 



(31), 



I O X -A- 1, 



pi^ovided the expansion should prove convergent. [Heine, §21.] 

A thorough discussion of the function H is given in the next chapter. 
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Chapter III. — On the Function H«(^). 
The fimctioB. H^^ (x) which, when mod. x is not less than 1, may be written 



H^^ [x) 



2 



u 



n(ra-2)n(-i) 

n (n - i) 



.^-« + lJ^(■!L_i,!^^^ + l ^-2) . . (28), 



\ 



2 



may be identified in any expression, in which mod. x may be as great as we please, as 

1 

that solution of (7), which vanishes when mod. x is infinite, and is equal to — ~ r: 

when X = 1. 

Thus, reverting to the expressions of p. 460, we see that, introducing a constant 
multiplier, 

n (tt - 2) n 



tl^ \X) 



n (n - i) 



^«-ip(_l,,,_l,^+l ^3^. . . (32). 



Again, by a well-known theorem 



11 (n — 2) I I (-1) -p / _ 1 



n 



li ^^ + h i^) 



w 







^l^ufil-^U^JdlL 



Hence 



H, {x) = -" t^' f t^^^-^ (1 - u)Hl «- t^f ¥ c^t^ 



• « . » 



(33), 



a form which may be compared with (11). 

In (32) and (33) no restriction is placed upon the value of x, and n is supposed to 
be a real positive quantity, but not necessarily an integer. 

From (33) we get a form analogous to (14). 

For write u^=. x — ^{x^ — 1 ) cosh ^ ; this gives, on substitution, 



H;, (itj) = — J d4 {x^ -. 1) sinh^ <^ {x — ^{x^ — 1) cosh (^}^^~^ . , (34), 



where 



. 11 ^ + 1 
9o = ilog^3I 



The sign which we attach to the logarithm is immaterial ; I shall take that sign 
which makes the real part, if any, positive ; or, as we may write it, 

9o = ¥lo8':; — ^, = log p •- ^a•, 
p > 1, and ^ i.'^ ^ ^- ^ 1 



^7T < O' < ^TT, 



iJ \J 4^ 
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In (34) write 



{x — ^{x^. — 1) cosh (f>) (x + y^{x^ — I) cosh ^) =: 1 . 



Substitute^ and we get 



XT /^\ -_ _ f 



CO 



(x^ — 1) siiili^ d dd 



(x + v^(?~^ 1 . cosh 0)^+1 • • • 



♦ » { tj %} f 



In this expression, x may not be negative and greater than unity. 
The expression (34) may be made to furnish us "with another, on the model of 
Mehleb's form for P^^ (x). 



For, write 






.^/{x^ — 1) cosh (J) := e^. 



f-^-^ 



Then 



^{x^ — 1) sinh^ = %/ie^^ ~^ '^xe^ + 1)? 



where the square root on the right has such a sign that its real part is of like sign 
with the real part of \/{x^ -«^ l ) sinh <^. 
Thus 



XT /..A 



cZX\/{2(coshx -- a;)}e("-*'^ • ■ • • (36), 



where 



'00 



sinh (cosh ^ a?) = — ^{x^ — 1 ). 



This expression corresponds to (15), p, 458; (15) and (36) may be differentiated 
once with respect to cr, but not twice. 

We have, by p. 465 , where n is a positive integer, 



n = oe 



"1 • "I -3 

t nin — l) (2n — I) H« (u) L (x) — -^^^^^^ + sc = — ~- 



« ft 



\.0 ( h 



HencOj multiplying by y 



■7* 



^ J-O iX) 

ft S CO 

t n {n — 1) 



n =: 2 



{2n — l)y lln((y)ln{x) -- II;, {y) {{n + 1 ) I,, + 1 (i») — (w --- 2) l,,^i {x)} 



since 



(2s — 1) 0?!,, (a?) = (w + l).I« + i {x) + {n — 2) %^^ {x) . . . (16). 



Hence we get 

(n + 1) H^j+i (y) — (?iW.>-l)^H;j (y) + ('^^ --2) H^^^i {y) =^ 



« (38) 



3H3 (y) — 32/H.3 (y) + 1 = . . . , . , . (38a), 
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A proof of (38) which, is valid when n is not a positiye integer ^ may be obtained 
from (36). 

From (38) it follows that 

where S^,, e^^, ^^^ are the quantities so called on p. 460 ; with the following exceptions 



xm^ (x) = S3H5 (cc) + egHg (x) — 1^1 



X-Hj (x) = 8^Ki {x) + e^Ha (x) 



If Q„ (x) be the Zonal Harmonic of the second kind 



th \ 



9 



• * » < 



Q,^^ (a?) = (l<^ {^ — y^{x^ "^ 1) <30sh (f>y^ dfj}^ 



where 



^j^ = ^ log 



^ +. 1 



"Heine, | 36.] 



(39 a), 



From this may be derived, as on p. 458 for Im{x)^ the expression 



rcosii-^A' 



where sinh (cosh ""^ x) 
Hence 





— 00 


-y{x''-l). 


dSn 


(^) 



^(»+i)*Vl% 



V^{2(cosh%-.4} 



^1^, - Q^^-i (^) • 



• t • » ♦ 



^it; 



. . (40)? 



and also 



{^Mi ~-~ i j xi^^ ^it^j — Ccj;^ (^j -— 'y^^^g (£iJ) ...... (41), 



Again, as on p, 461, we find 



dXl vl 



Xl;^ \XJ "^ \:^u \^) "T -^^M^ — 1 \*^/ "^ ^}?*— 3 («^/ . • . » . (42). 



HI 'Y* I * " """ " *' 



nov- 2)n(- L) 



2«n(«,-i) 



'■ Di(«) (a;) , 



where Di"^") (a?) is defined by Heine, § 3L 
By means of the sequence equations 

(n + 1 ) I«+i (^) — (2w — 1) xl,, (x) + (fi. — 2) I^,_j (a?) = . . . (16) 



(n + 1) H/, + 1 (a?) — (2ti — 1) ^H^^ (^) + (^^ — 2) H^i^i (x) = . . (38) 
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we see that 






> 



* • 



(43), 



.J 



where ar, fir ai*e rational functions of degrees 7\ r — 1 respectively in x. This holds 

whether n is integral or not. If 7i is a positive integer, we have, taking I'egard to 

the equations 

3I3 (x) — Sxlc^ (x) =0 

3H3 (x) - SxRq (a?) + I = 



K„{x) — a^.gHa (a?) + /3,,„2 



e » « 



• » 



But by (34) 



Jo 



• I ^ij A. It 



X 



|l0g(.^^+l)/(2;-l) ,y. ^ I 

(oj^ — 1) sinh^ ^.d^= -| I3 {x) log^^^J — - — j . (44a). 



Hence we get the expression, n being a positive integer, 



ri;^ y^} — 2 ^u \p^) ^^^ 



dz — - 1 






* i t » 



(44), 



where K^^ (oj) is a rational integral function of x^ of degree n -— 1. 
By substituting (44) in the equation (7) 



(1 _ x^) ^# + ;^ (v^ ^ 1) H, (a^) - 



dw 



we can obtain an expression for 'Kn{x) in terms of the l^s, viz. 



\ / (^2r — 1) {n — 7') 



1 



(2r — 1) (^i — r)' 
n {n — 1) 



l;^_2r+i • » {44b), 



the series beginning with Iq or Ij. [Heine, § 26.] 

The expression (44), where the logarithm is defined on p. 467, might be taken as 
the definition of H,, {x) when n is a positive integer, without restriction upon x ; but 
when X is real and less than unity, it is advisable to modify it ; viz., from (44) we get 



1 1„ {x) log YZ", - K» (iK) ± V • I» i^) • 
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In this case, I take 



1 _L. QQ 

H,, {x) = 1 1,, {x) log Y~-~ — Kn (a?) ..... . (44c), 



m 



a form which^ along with ln{^), is capable of expressing any solution of (7). 
The form (44) is easily derived from 



H,, (.0?) = I J^ -"^^^ ....... (44d), 



■1 ^'^ - y 

which itself follows directly from the development (37). [Heine, | 28.] 
The two forms (SS), (11), are simply transformations of the two solutions 

y = [ 7/-1 (1 - i^)v-^-i (1 «. tu)-'' dth, /3, y — A positive 
qj = Y li^"^ (1 — 1^)^""^"^^ (1 — . tuy^- du^ y — ^, 1 — a, positive 
which, among others, Jacobi gives (' Crelle/ vol. 56) as solutions of 

t{i - ty^. + l7 - {<^-h ^ + i)t^% - o^^y = 0. 

See (7d) p. 460 ; m = t^ — 1 in the first, and = — ^^ in the second, t = £^, and Jacobi 
supposes it positive, but here that is not necessary. 

Let us now coordinate the various solutions that have been obtained. 

The different series used are only permissible expressions for the solutions of (7) so 
long as they are convergent. 

Now F (a, /3^ y, x) is convergent, 

when <% -j- ^ -^ y IB not positive, if mod. x < 1 
when a + j8 •— 7 is negative, if mod. x ^ I. 
Hence 

(9) (28) are permissible for mod. x '^ 1, 
(9a) (9b) are permissible for mod. x 'jf> 1, 

(10) (32) whatever x may be. 

Now, by direct development, we know, (^, y — . ^^ positive), 

I V~^ (1 ^ u)y-^'-' (1 -^ xu)- du = '^^^:zj}I^ p (^^ 0^ y^ .^y 

Hence (33) = (32), and, therefore, = (31), when this last is a permissible expression. 



472 



MR. E. A. SAMPSON ON STOKES'S CURRENT FUNGTION. 



But (11) ^ (9), or (10), unless wis a positive integer. For (ll)s(l4) and from 



(14) 



Xjv " 



J-» \p) 



x~l^ 



^A 






(9) 



^" — 1 



2 



n,-l 



n (fi >- f ) 
n (^i) n ( "^^4) 






«, 



1 



9 > 'i 



r?., X 



2 



(10) 



Ijf^ \^/ 



■I 



S 



n (7t - 1) 
n 0^) n ( - i) 



' M + J 



-F( 



9, 



^^ 



3 3 

? 2? 2 



". f )■ 



But the series on the right are divergent for x = L Similarly (9a) and (9b) do not 
agree with (11). 

Taking the definitions (11) and (33) for I,^ (^) and H,,(o?), we find 



I,: (x) == e''"" I,, ( ~ x) ■ H,, {x) = e/ ^'' ^ ^> B, ■ ( — x) 



» » '' » c 



(A) 



I4±l) = H«(±l) 



n {n — 1) ' 



» « S 3 



(B) 



In (0) 



2^(J?^--|)^(i) ntr 

TT n(-|W') * ' 2 



H, (0) = ^^^^^-^^^^^^^ ( , cos 



sm 



nir 






(C) 



.lib t!C X.JI \X) — - 2i 



« — 1 



fCss CO 



n (7t -- 1) 
nwn(-i) 



JLju «X" .JlI.^^, \*^ ) 



:c = 00 






(D) 



Add to thesOj when ^^ is a positive integer^ 



n odd 



n even 



a; — '^^ ■ 



P_, (0) = (-!)'(«- 1) 



n (I «. - 1) 



»-]. 



n(i«--i)n(~i) 



Q„_i(o)=(-i) 



. , n (i ^t - 1 ) II (i) 

11 a n - i) 



)■ 



J 



» • 



(E) 



EINE, jil 4, 25 J 



These results will be quoted as , , . . 



a A 9 i » « ji 



• (45) 



We can now give the expression for H^^ {x) in ascending powers of x. For 
the expressions (9a) and (9b) satisfy the equation (7), whether n is odd or even ; and 
where I,^ {x) can be developed in even powers of x, H« {x) can be developed in 
odd powers^ and conversely. The constant factor is detennined by (45) ; hence 
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n even 

.1.1 \ ^'v "~~ *<> / \ lU Zi I 



n odd 



H,, (a?) — (— 1 )*('^* + i) • ^ ^^ ^ f* / -. ^^ ^- ^ 1^ 1^ ^a j ^ ^ (28b), 



These apply when mod. x\ \\ and only when n is integral ; for we see from (45a) 
that neither I^^ {x) nor H^ (x) can be expanded in a series of integral powers of 
X when n is not an integer. 

Two functions, R.^, T,,,, have been used by Mr. Hicks^'' in connection with his 
researches on the motion of circular vortices, which are such that 



R^ = — -- — I^^_^^ (cosh te), 



4^3 _ 1 
n — A ijL^ .{.i \^COSn %)* 



Mr. BASSETt employs a function L^, where 

2 

and finds the sequence equation corresponding to (38), and the values of Lq, L^ in 
elliptic integrals. 

Mr. Hicks gives equations which are equivalent to (18) and (40); approximate 
expressions for the first five R's and first four T's in terms of ^ = e~^ (or, if we 
write X = cosh u, k is the f of p. 456) ; and an expression for T^, which leads to 

fir 

Iln + i{x)=\ Gosn0^{2{x~ Goa0)}de (46). 

A proof of (46) will be given in the next chapter (p. 483) ; but I will here show that 
it applies only when nis a. positive integer. 
For 

["IT 

H^+i (x) + H^^i {x) = 2x11^^^ (x) — 2 v/2 cos nO {x — cos Of dd. 

* 'PML Trans.,' 1884, 1885. 

t * Hydrodynamics,' vol. 1, § 115. 
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Now 



2^/2 cos n0 {x — cos 0f d0 



2\/ 



9. 



n 



sin nd {x — cos &f 



w 







+ 2^{H„^,,(a;) 



-O,??— .-I \*^/ j 



NoWj ^ B i^ a positive integer^ the expression in 
we find the sequence equation 



] vanishes at both limits, and 



{n + I) H,,+| (x) — 2nxll,,+^ (x) + (w — |).H,,^j (x) = , . • (88). 
We see from p. 470 that we may write 



-fci^^^i {xj = (iff,..„ikl^\X) -~" Pu^ixlx \Xfp 

where^ n being a positive integer, et^^-u A»-i ^^^^ rational integral functions of a? of 
degrees n — 1, n '— 2, respectively ; and L, L, Hs, Hi may be expressed in terms of 
complete elliptic integrals of the first and second kinds. 
For we have, by (14a) 



1| {x) — 



and 



where 



1 r^ 

TT Jo 



IT 



(oc^ — 1) sin^ dd 
(m + \/a}^ -— 1 . cos dy 

fir y^^,2 _ I ^ QQg g0 




{m + \/# — 1 COS #) 



Jo 



cos 



6 \/x 



B 



1)* 



1 f 

^ 4- cos ^ ^{x^ •— 1) = i 1 



^' 



3 



:s 



ti/ 



{m + oos0\/^^ — !)*_ 



>d9^ 






% ^ 



2 



Also, by (14) 



and 



It % 



2 

TT 



^.E(i'/,)-,.f•F(F.f 



i| (i») 



1 ■ fir 



{gc? — l)siii^#cl# 
(a? + cos ^ V^^^ ■— 1) ■ 



s ♦ 



2 P 

w J 



V^i^""^^ cos # (a; + cos 9 's/x^ — 1)* dd^ 



{x^^l) sin^ ^ =2 -- v'aJ^ — 1 cos ^ (i» + V^^ — 1 . cos 0) 

-^ x{x •{- %/ x^ — 1 * cos ^) — 1 



, (47). 
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or 



m fir 

•. I| {x) = — |I.| {oc) + - (x + COS Q V a?3 — 1)* M 

TT Jo 



1 



fesym 



TT 



dd 



{x + cos d -^ x^ — 1) 



-'-« w 






37r 



x^-^mi'. 



|M»Mr 



^'F A;', 



^ 



r 5 



• * • 



(47a). 



Also, using the results given by Mr. Basset ('Hydrodynamics/ vol 1, p. 108), we 
have 



Hi (jr) 



u 



t* E ( ^, 



TT 



IV(«='-1)FK' 



TT 

2 



H, {x) = -i [r* 0= E (i, I) - r* v'{x^ - 1) F (^^, 



TT 



* * * 



(48). 



The expressions (48) give very little new information about the functions H, since 
we have already obtained in (32) the expression for H^^ (x) in an ascending series of 
powers of £; but (47) throws some light on the difference between the forms (10) 
and (11); for having regard to the expressions (Caylby's * Elliptic Functions/ 
Chapter IIL, 177), 



TT 



F| y/ 1 — f ' 5 



4 12 
log f + ^, f 



'2 



22 



1JC 



1.2 



12 92 

1^ 09 A9. ^ 



2^.4^ 



log 



i.2 



3.4 



-^ &c.. 



w 



E( v/1 - e, ~ 



1 +i$^ 



1 Tit JL 



JL « iM 



1^ ^ 

™ 22.4^ 



4. 

02r T. 



2 






— 1 
3.4J' 



it is at once obvious that (U), which leads to (47) and (47a), gives for I;^ (x) a function 
of £ which, in general, is not identical with (10). 

We shall have to deal with infinite series into which I's and H^s enter ; it is, there- 
fore, requisite to know to what values these approximate when their order, n^ becomes 
infinite. We have, whether n is integral or fractional, 



TT (nr) 



TT /m ^ X\ ^ \^^ ^^ ""' ^ ^^ » ^' ^ / • • • (^2). 



whence 



11 ('?!■ -g) 



-tl/* \«^/ — 



IT' 






• • » • • iTtt/ f , 



when n is very great. [Heine, § 40.] 
Again, 

t—n rl 
1^ (x) = - - t^'^^ (1 — t^)* (tt 

3 p 2 



1^)^ (ii^ 



t • 



» • ♦ » \i-j-./« 
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Now, when n is very great, each element tf ( 1 — uf (u — ^^)* du of the integral is 
very smallj except when u is nearly equal to 1 ; we may, therefore, pnt 



I^ (x) = ^ — !^ — -^-^ tif^ ( 1 — uf du 

TT Jo 



when 11 is very greats 



"^ TT n 0^ -t- 1) ^ ^ ^ ' "^ 2v/r7r#)^ ^ ^ ^ ... \^yA;. 

The expressions (49) and (41) a) enable us to discuss the convergency of the develop- 
meno 

- — " — - = % n{n -— 1) (2n — 1) 11,, [y) 1^ {x) (37). 

For, from the equations (16) and (38), we find 
%n{n --- I) (2n — 1) (y — x) 1,^ (x) H;^ (y) 

2 

= 1 -- x^ + {n -- 1) n {n + 1) [L(a^) Hn+i (y) — H,, (y) ln^.i{xj 



Now the limit when ti is very great of the expression in [ 



IB 



4Mi)"^<'-f'-'-''')-('-i. 

Hence (37) is true so long as mod. rj > mod, ^. [Heine, | 45.] 



Chapter IV.— On Functions appkopriate to the Spheroid and the Tore, 



We have seen in the introduction that the expression 

tfj = (Ar^^ + Br "-« + !) {CI., (cos 0} + D H^cos 0)} 
where A, B, C^ D, and n are arbitrary, is a solution of 

and we shall find that sokitions expressed in this mam i or are appropriate to the 
discussion of questions of fluid motion, when the boundary is a sphere. 
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In the present chapter I propose to develop solutions which shall be appropriate to 
the spheroid and the tore. 

Consider a spheroid whose centre is at the origin, and whose major axis, and axis 
of symmetry coincide in the axis of z. Let the distance between the foci on the axis 
of z be 2h, Then if ph, qh, be the major axes of the confocals to this spheroid which 
pass through the point (xat, z), p^, q^ are the roots of the equation in X, 



.3 



So that 



K — i A. 



z =2 hpq . 
OT = ih Y^( 1 — 1>^ . 1 — q^) 



(50), 



where we shall suppose p^ algebraically greater than q^, and, by an arbitrary con- 
vention, q to change sign with z. 
From (50) we find 



dw^' + d^ ^hr (q^ -— j^^^ 



dp^ 



dcf 



1 



p 



f 



whence the normal distance between the two spheroids jp? i> + dp, is 



^ V (iTrli) ^-P' 



and between the hyperboloids q^ q_'\' dq^ 



h \/ ip^-^^ ) dq. 




Now, reverting to the hydrodynamical origin of the subject, it is shown in the intro- 
duction that — tjt""^ Di// is the circulation per unit area for any closed curve drawn in 
a meridional plane. Consider the small curve bounded by p, p + #, g, q + dq. 




The velocities in the directions of the outward-drawn normals are 



478 



MR. R. A. SAMPSON OK STOKES'S CURRENT FUNCTION. 



to p, 



P = 



1 /I — q^ d^lr 

im V ^^ — • q^ dq 



to q^ 



Q 



h 



1 /I — p' dyfr 

;.nr V q^ — p^ dp 



I 

> 



(51), 



-J 



and the circulation in the circuit considered is 



ttMJaaa 






(ip (ig'j 






But this is equal to 






c^2^° 



1 — ^'^ dp'^ 1 -— ^^ {^g'^ 



rZp dq 






hence 



Di// 



A^ (■^^ — 2'^) 



and the equation 



transforms into 



{i-f) 



Di|;=0 






(1 - q^) 



d^y]/^ 
d(f' _ 



• « 



e e e 9 » 6 » 



V • 






. 



« o « « e e 



(6a), 



(3) 



(36). 



In obtaining the above I have supposed the ellipsoids and hyperboloids of revolu- 
tion about their transverse axes ; and, in this case, p, q are the reciprocals of the 
eccentricities of the generating curves. Hence, for ovary spheroid, 

p lies between 1 on axis, oo at infinity ; 

q ,, 1 ,, , at plane of sj^^mraetry. 

If we wish our equations to apply to planetary spheroids and hyperboloids of one 
sheet, we must suppose — h^ to be the squa.re of the radius of the focal circle in the 
plane of symmetry, and p, q are equal to (e^ — 1)*M where e is the eccentricity of 
the generating ellipse or hyperbola. 

Hence, for planetary spheroid, 

p lies between within ring of foci, ico at infinity, being everywhere a pure 

imaginary, 
q lies between on plane of symmetry beyond focal circle, 1 on axis. 



Hence, whether we divide up space by means of ovary spheroids and hyperboloids of 
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two sheets, or by means of planetary spheroids and hyperboloids of one sheet, q is 
real and lies between and 1. 

Hence, by p. 462, any function of the coordinates, which remains everywhere finite. 

may be expressed by the series S In {q)fn (p) where n is a positive integer. If this 

% = 

is a solution of (3) we find on substituting 



% = 00 

S In (q) . 
whence. 



(l-f)'^ + n{n-l)fAp) 



= 0, 



m^ 



which is the equation (7). 

Hence, any solution of (3) may be written 

t In {q) [A J^ (p) + BMn (p)] 

where ^ is a positive integer. 

When ti is a positive integer I^ (p) I^ (q) is a rational integral function of the 
coordinates t^r, z, of degree n. 

In the first place, let n be even. The product may be resolved into a series of 
products, {p^ - a?) {q^ - a^) where a is any root of I, (x) = 0. 

And by (50), 

when a = 1 ; 

when OL zjz \, 

But by (9) we see that, if n (a^) denotes the product of all the a^B, 

n (a^) = - ~^'"^ , 

^ ^ 2.4...7^.7i — l.Ti + 1...2^ — 3 

and 

n (1 -- a^) = Lt M± ^ ElzjSH = L 

Hence, when n is even, 

l^pjl^g)- 2.4....7..^~1...2^-3 2A^-^H^'^l-"«' ^j. • i^)^;. 

where all the roots occur in the product excepting a^ = 1. 
Similarly, when n is odd. 



T«(P) ln{q) - 2A..M -^ l.n.n + 2,..2n ^^^ m^^^^Xa^'^ 1 - a^^ ^^ 



where, under the product, all the roots are included, excepting a=: ± 1 and ct =:^ 0. 
the latter of which gives the term pg, or zjh. 

Now it is easy to see that, except for a constant factor, there Is only one homo 
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geneous rational function of the coordinat6)S of positive degree r which satisfies (3). 
For such a function in its most general form involves r + 1 constants, and, when we 

have operated upon it with ^{^r ( -—^ + Ti "^ T ' ^^'^ ^-^^^^ ^^^^ with a similar function of 

degree 7^—1, each of whose terms must vanish. This gives r linear equations 
between the r + 1 constants, showing that one, and one ouly, is independent, and 
that appears as an arbitrary factor. 

Hence the r functions I^. (p) I^ (g^), '^r^i{p)\^i{([), - • • must be equivalent to these 
r homogeneous solutions of (3), and are therefore sufficient for expressing any rational 
integral function of the coordinates of positive degree r, which satisfies (3). 

We shall require hereafter the solution of 

where y^(_p), <^,, (g) are linear functions of 1m, {p) a^nd H.,n{p), I^, ($') and H^(g'), and 
m ::^ n. 

A particular integral is 



^h 



_- A (y) ^n (q) 

(n — m) (n -{- m '-- 1) 



I shall now consider the functions which are appropriate to the tore. 

Write 



trr + ^2; :=: c coth 



2~" 



Then 



sinh 7] -{- i sin ^ 
cosh 7} — cos ^ 



. (53). 



0, 



^^ ^ z^ ^ 2zc cot 1^ — c^ = 
zs^ + z^ •— 2otc coth >7 + c^ == 0. 




And if P be any point in a meridional plane, A, B two points in that plane sy 
metrically placed at distance c on opposite sides of the axis of symmetry, 



m 
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eV^ 



f 



BP 
AP' 

APB. 



Hence the surfaces represented by 77 =:::: const, are tores, and those represented by 

^ ::= const., the surfaces given by the intersection of two equal spheres. 

We also find 

dTf + dp 



d'ns^ + dz^ == c 



3 



(cosh 7] — 00,3 ^y 



Hence the normal distance between the surfaces 77, 77 + di] is 



and between f , f + d^, 



6' ^^'/; 
cosh -t; — " cos f 

^^^^ 

cosh 7] — cos f 



ThVis the velocities in the directions of the outward-drawn normals to the surfaces 
rj, i are 

cosh 1] — cos ^ chfr 
cm 



H 






cosh 7; — cos 1^ cl's^ 
cm drj 



> 



.y 



and hence the circulation in the circuit formed by >?, -/? + d-q, ^, ^ -\- d^, 



(54), 




fL 
d^ 



H 



cosli 7] — cos I 



(i^d77 + 



d 
d7] 



/— ^ 



and we get 



where 



cosh 7j ^ cos ^ 

d "^1 d-xfr 

d^ [m d^ 



+ 



d^dr], 

d 



d7] 



1 ^:t 

m drj 



Bxfj 



(cosh 7] — cos f )^ 



7aT 






I«T 



c sinh 7} 



MDGCCXCI. — A. 



cosh 77 — cos ^ 
3 Q 



d^dyj, 



{6b) 
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If we write 



%jj = ^/(cosb. 7] -— cos i)^ 



<s 9 & m ^ 



« » * ItJOl* 



and we find 



1 d'^ 

1 ^l|f 



Bx}^ 



c sinh '?^ 

1 

c sinh '?; 



(cosh ?? — cos i) 

9 

(cosh 7) — COS ^) 



1^ _ % ^^^ f 

c^l 2 (cosh ?7 — cos I)* 
r% 2 (cosh •?; -— cos ^)*^ 



"^ 



V 



mn^ 



(56), 



(cosh^'?7 —'"cos f)* 



C' 






coth ^ $ + 1 



7} 



4 



(^ — COS ^ 






(Z' 



.^r 



(1 - *'^) X.2 + i 



CvJb 



dri 

{x — cos I)* . t^ 



, (6c), 



where 



and 



X 



cosl 



iX f / « 



I)i|»= 



© * 



««« tt««»# 



(3) 



transforms into 



where 



^ /I ^2\^ I % 



■^:=: \p (cosh 7J — COS ^)' 







» • • » 



(do)j 



Now ^5 as it occurs in the equation (Sc), is a real angle lying between w and 0. 

Hence, any function^ x> which satisfies (Sc)^ may be expressed^ by Foueiib's theoremj 

in a series 

tfn (x) cos (ni + a,^), 

where 7i is a positive integer. 

Substitute in (3c), and we find 



If — cos (ni + oCff) 



w^hence 



dj f ('iF) 



0, 



which shows that 



(1 - .■^^) ^-^~ + (« + i) (» - -i)y« («) 



3 « • » 



(7), 



B being a positive integer, and any solution of (3) may be written 



t|^ 



{x — cos f )* 



2 COS {ni + «"#) [A,,I;,+| (a?) + B,^H#+i (^)] . . . (57), 



an expression identical, except as regards notation, with that given by Mr. Hicks. 
(^Phil Trans.,' 1884.) 

We are now in a position to show how the expression given on p. 4^3 for H;,^|(aj) 
may be derived, [(7f, Hicks, 'Phil Trans.,' 1884.] 
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For ^|r = \/{^^ + ^^) is a solution of (3) ; and 

cosh 7} — cos ^ 
Hence we may write 



n := 00 



{x + COB if .:= t COS {n^ + a,,) [A J,, + . (^) + B,,H,, + . (x)]. 



n-~ 



But it is clear that the left-hand member may be developed in cosines of multiples 
of f, and, therefore, a,^ = 0. 
Hence 

A,J,,+x (x) + B;,H,, + x (x) z=z -^ COS n^ ^{x + cosf) (if. 
Now, when x =■ co ^ the right-hand member vanishes. In fact, since 



TT 



cos nf cos"^ ^d^ = 0, 





when n, r are positive integers and r < n, it is clear that the highest power occurring 
on the right when the expression is developed in descending powers of x, is x'^'''^^. 
Hence An= 0. To find B,„ put x = 1. 



But 



f^cos ni y2cos ldi= ^^"J^(;rr^^ 



H^^^ + i(l) == — 



therefore 






B,, — 

TT 



and 



H.+^ (^) = (- 1)- r cos n^s/i^ {x + cos f)} c^f 

Jo 

= G0Bni^{2(x — COS i)} d^ (46), 

where n is Sb positive integer, and x is real and > 1. 

These functions have been employed by Mr. Hicks in the papers already referred 
to ; and, except for purely mathematical interest, it is chiefly in connection with the 

theory of fluid-motion about a tore that functions of fractional order claim our notice. 

2 
Mr. Basset'''''' has, however, employed the function L,^ == -■ (— 1)^ e*"*^^^ H^ + ^ (cosh rj), 

to find the current function for motion in an infinite liquid due to the rotation of an 
infinite cylinder whose cross-section is a lemniscate. 

* ^ Hydrodynamics,' vol. 1, § 115. 

<D ^4 -*-«/ 
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Any solution of (3), which consists exchisively of positive integral powers of the co- 
ordinates OT, z^ can be written in the form 

1 '/l ;= CO 

S B;,H,, + x {x) cos {n^ + a,,). 



To arrive at this theorem we notice that 



But 



hence we see that 



7 

2 i'^ ix — cos ^\^' -^ z (x — cosf )^\ 
2 -;-! — -™ -sy^ v^ 4- (^'^ "- ^'^ "^ 1) ~r ; 
2M ~ j (o^. — cos ^f ^ U,. (a? — cos f)% 



where U,. is a rational integral function of z^, z, of degree r, and consecutive functions 
U^,, JJr + i are connected by the relation 



„7TT ^TT 

U,^.^^ ™ ~- 2wz "T-^ + (^^ -- z^ — I) --jf + 2:11,, 



3 



and Uq = ] . 

Forming the first five fanctions, we have 

U] ^=: Z^ 

U2 :^ ts-^ — Ij 

U5 :^ z (45^^ — 60?^%^ — 30ot^ — 1). 

From these it is clear that when r is odd, 2; is a factor of U,. ; but that, otherwise, 
U;. is an even function of the coordinates, zn.z^ — m occurring in higher power than z. 
Now 

n=oo 1 

(x — cos f )*, — t - H^^ + i- (^) cos n^. 



n 



= 'TT 



Therefore^ 



(,^' — COS f> ,,==0 'TT -^^ V / y 2 / 



and therefore xfi = JJr is a solution of (3). But there are only r rational functions of 
the coordinates of positive integral degree not greater than r, which satisfy (3), and 
which are independent. Hence, in terms of the functions, JJr, U^'-i, &c.^ we can 
express any such solution of (3) ; whence we obtain the theorem stated at the 
beginning of this section, and we see that B,^ is a rational integral function of n. of 
the same degree in n as the solution of (3) is in the coordinates, 1^, z. 
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Table I. — ^Functions of w and z of Positive Degree satisfying 



dm'^ dz^ m dm 



(1) Z, 



w » 



(2) 

(3) zm^ 

(4) OT^ (42^ -- 7;y2). 

(5) 2ot^ (42^ — 3zs^). 

(6) ^3 (82* — 12^^ + S7*). 

(7) Z7s^{8z'^ — 202^^2 + 5t5*). 

(8) ts^ (642" — 24O2V + 120zV — 5^"). 

(9) 2w«(6428 — 336zW + 2802^^* — 35m'^). 

(10) z!y3 ^12828 - 8962V + 11202*w* - 28O2W 4- 7zs% 
The above were calculated from the formula 

U,, = {2n ■— 3) 2U„_i — (ar^ + 2^) "f~^ ' 
which may be readily deduced from (13). 

TjiBLE II. — 1„ {x) in powers of x. 

w = 0. — 1. 

rvti , -^ I 



3. 



A 



(L/ • 



6. 



2 


■0 




^3 — 


0? 

* 




2 




5.'^* - 


- 6a?2 + 1 






8 




7a?^ - 


- lOx^ + 3a; 






8 




21^-6 


- 35«* + 15a;3 - 


-1 




1.6 




330^7 


- 63a;5 + i^«? - 


- 5x 



8 



9. 



10. 



16 

429a;8 - 924a!'' + 630ai^ - 1403;^ + 5 

715a;« - 17 16a;7 + 1386iB^ - 42 0a;^ + 35a; 

128 

2431^1° - 6435a;8 + 6006a;'' - 2310a:^ + Slda;^ - 7 

256 
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Table III.~K,^ (x) in powers of x. 



n= 2 



X 

2 



3 



^ • I ^4 _ 2.3: o^3 l_ _2_, 

6* -^JL /y5 _/_ /yiO I. Jt-JLjti ^» 

' 16"^^ 4 "^ "i 2 4 "^^ 

3 3 ^6 1 3 ^4 L 1 3 ^2 8 



105 



»7 . 3. 3. ^6 JL 3 ^4 I 1.0 3. ^5 

/ ? 16*^ ""■ 4 «^ "I 80 *^ 

Q . 4.2 9. ^7 7 8.1 ^5 J 2.0 3.£ ^3 1873 n., 

*^^ 128'^ 128*^ I 64 '^ "^ 4480 '^• 

Q . 715 ^8 4 4 3.3 ^6 I 9.5.2 a.4 _ 6.9 6.2 nr.3 I. _16_ 

*^5 128'^ "^^ 384"^ "Tl28»^ 4480*^ ■T"315* 

1 n • 2.4 S^i ^9 _ 8. 4 3 7 ^7 r %.1_\1. ^5 _ 4 3 6 7 ^3 I 3 5 6 3 ^, 

LU^ 256*^ 384^"T" 160*^ 896'^"r80460'^- 

These were calculated from the sequence equations (38), (38a). 

7^K,^ {x) = (27^ — 3) xK.n^i {x) »-» (^ ~ 3) K^._2 (r^), 
3K3 (ir) =: 3ccKp^ (iu) — le 

Table IV.— -ii;^^ in terms of I (.'r)^s. 

ti = ; — Iq, 

3 ; 2I3 + Ii* 

A- 8T»Ll2T_T 

5 ^ y Tg + -7- T3 + Ip 

" ^ 2 1 -^6 » 3 H "T~ 7 -^2 -^0^ 

^ 5 ¥3 I7 + TT I5 + ~3~ I3 + ^l- 

8- 128Tj_64T|112T|8 1 T 

^ ^ 4 2 9 -"-8 "T" 3^9 J^6 "T 33" ^4 T -3 ^2 "™' -^0* 

q . 1. 2. 8. T l6^4.T I432T |40T I.T 

^' 715-»-9"T"55%"r T4 3 ^5 + XT ^-3 + ^-1' 

1 • -2,5.6_ T ll.9 20Tr32Ti560Ti3.OT T 

^^ y 243 1 J^IO + 24 3 1 J-8 + 13 ^"6 + T43 H + iT -^2 H' 

Table V.~Cos (7^ cos""^ x) in terms of I (x)^^, 

n = ; — Iq. 

1 ' T 

AJ ^ 3 ~°"~ 0' 

3 ; 8I3 + Ip 

A ' J5 4 T I 1.6. r T 

^ ^ 5 J-4 "1- 5 -^2 """ -^0' 

n . 12 8. T _L 4^ T I T 

^ p 7 -'^S I 7 -*-3 n ^1» 

A . 5JL 2 T I 1. 2. 8. T 1 J^O 8 T ™. T 

^ ^ 2 1 ^6 r 1 5 ^4 r 3 5 ^2 -*-0' 

7 . 10 2. 4 T r X 2. 8. T r i 6. T 1 T 

' J 3 3 J-7 i" 11 -L5 "T 3 J^3 n^ ^P 

Q. 16 384 T l_ 4 0^ 9. 6. T I. 2_ 5. 6. T I. 6.4 T T 

°? 4 29 -*-8 "I" 2 73 -*^6'T' 3 3 -*'4 "t" 2 1 -^ -"-O' 

Q. 3.2.768T _|1024T _i 10 3.6 8T J_400T i_T 

^ ^ 7 15 H "T 5 5 ~ J^7 "T" 1 b OT" J-5 -f -7 7 X3 -f -Li- 

in. 1310 72 T _! 16.3.8.4 T _l 5.12,T 1. 3 2.0.0 f l lOOT T 

^^ ^ 2 43 1 -*-10 "T" 7 £5 9 3 ^8"!" 3 9 ^ 6 "I 4 2 <) H\ 3 3^2 % 
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Table VI. — Various Forms of I^ (x). 



on 



TT / 3 \ 

n(^)n(-i) 



fi % 



^^^l + h-^^ + h^ ^) • 



. . (9) 



(-. iy^/3 



n 

11 1 ^ — -^ 



F - 



-^l^ 



2n 1 ~ ) n ( 



2/ 



ft fl 



1 J /y»(i 



• • > 



(9A) 



(- 1)^^^* 



1) 



11 ( ^ - 1 



7h 



xF 



n(2~i)n( 



2/ 



-4-1 - la;^' 

2 ^'2' 



* « 



• (9B) 



n (Ti, - 1) 
nwn(-i)- 



cos {n cos ^ a;) — ;: 



li'/i 



2 (2n - 3) 



COS (n — 2 cos"-^ a?) — . 



. (lOA.) 



nJjJIllL ^-» F ( - 1 - ri, - % + f , 1^2) .. . (10) 



- u'^-^ il - uf {u - ^f du (11) 



TT J|3 



1 ^^-Y V^^^V ' (12) 



\n — 1 \dwj 



2 



/ — 1 y^ "" 1 ~ - 



(13) 



1 

TT 



TT 



s 







in^ '?; {x^ — 1) (^ — ■ cos ')y s/x^ — 1)^* ^ d^'jy . 



(14) 



- I sin^ 4> (x^ —l){x+ cos (^ y«® - l)-»-i (i,^ . . (14a) 



1 rcosu "** '' ft/' 

- c?.^ y{2(x- cosh <^)}e (''-*)* . . . . (15) 

TT J — cosIi~"i X 



cZ^ •v/{2 (cos <^ — a?)} COS (^ — :|) ^ . . . (15b) 

TT J 



d(l> y/{2 (x — cos (f>)} sin (^ — -|) ^ . . . (15o) 

J cos — !« 



2 T'r 

'TT w cos— *^ 
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Table VII. — Various Forms of H„ (x). 



2 



_,^n(^-2)n(-i) 

n (n - 1) 



, , _, /m, — 1 n 



1 ^ ~" 2 

2 5 



* « 



(28) 



11 



{— 1)"'' 



tt/'Z 



n ( - - 1 ) 11 (i) 



11 



n 

2 



2 






e o 



(28a) 



(-1) 



U^^tl) 



n 

9. 



3 

'2 



n (i) 



fTi 



n 



n 



■11 n 



\ 



1 1 rv.^ 



9 • 



(28b) 



lUn - 2) ri j|) 



f-^F(-i,n-~-l,n + i,a . . . (32) 



n 



-^ 



n " ] 



Vj 



n 



■2 







(1 — uf (i — 'i^P)* (it^ 



« 



e » 



[ 



log {x + l)/(a; ~ 1) 



(33) 



d<l) (a-^—l) sinh^ ^ (a; — -v/x" — 1 cosh ^)"-2 . (34) 



\ 



CO 



- j £^6* (sc^ - 1) siiih^ ^ (« + ^/stf—l cosh ^)-«-i . . . (35) 



'cosh — i?j 



c^^ y{2 (cosh ^-ic)} e («-*)* (36) 



— 00 



x-\~l 



r I \ {x) log ^ _ ^ — K,, (iX') . 



e 9 « * » 



At. 



i I.. (^) loj 



1+0^ 



Ja.j,j{ \X 



J. 



x~ y 



r + 1 



e s 



ace 



v_ 



. . (44) 



• > • • ItcjCV/'/ 



, . . (44d) 



rcos(ti-~.|-)6> y{2(x — cos6^)} d^ (46) 

Jo 



(9), (9a), (9b), (10), (10a), (12), (13), (15c), (28a), (28.b), (44), (44o), and (44d) 
require n to be positive integer. (46) requires n of form integer + -|. (9a) and 
(28a) n even; (9b) and (28b) n odd. (28a), (28b), mod. x :|> 1. (15c) and (44c) 
X real and < 1. (28) mod. x ^ 1. (46) x real and > 1. (14a) x not a pure 
imaginary. (35) x not real, negative, and > 1. Subject to these limitations, 
the above forms give the same values, with the exception of (44c) among the 
functions H. 
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Chapter V. — Irrotational Motion in Hydrodynamics. 

The chief interest attaching to Stokes's Current Function lies in its application 
where the motion is rotational, and, consequently, where the Velocity Potential fails us. 
But it is never without value, since the equation xp = const, gives the stream surfaces. 
Moreover, methods are best illustrated by their simpler applications, and the results 
serve for comparison with more complicated cases. Accordingly, in this chapter I 
find the value of \}j for certain known cases of liquid motion. 

The problems we can attack with Stokes's Current Function are those of motion 
symmetrical with respect to an axis, and this requires that the boundaries shall be 
symmetrical. The boundaries may, therefore, be such surfaces as a sphere, a tore, a 
hyperboloid of revolution of one sheet. We may suppose the boundary and the 
liquid at infinity, either at rest or in motion, symmetrically. When the liquid at 
infinity is in motion, I shall suppose its motion uniform ; and we may then reduce any 
of the cases within this range to the problem of flux past a fixed obstacle, of liquid 
moving uniformly at infinity. 

The motion being irrotational, we have D\j) = 0, and consequently \jj consists of a . 
series of terms such as we have found as solutions of (3) in the foregoing chapters. 
If these are expressed as functions of the coordinates z^r, z, it is easy to see that the 
only positive power of the coordinates that can occur is zr^ ; for any higher power 
would make the velocity infinite at infinity, and the term z would make it infinite at 
the axis. This immediately strikes out various terms. For example, 

r cos 6 ; r^ I3 (cos 6) ; r^ l4(cos 0), &c.. and I;^ (p) I^ (q), I3 (p) I3 (9) . . . by p. 479. 

Our process is, then, to assume a value for \p which contains a complete series of 
terms, omitting these impossible ones. The mode of expression that is appropriate 
differs with the form of boundary. We then use the principle that x]j = const, is a 
stream surface, and may, therefore, be made to represent the surface of the obstacle ; 
while at infinity i// = — • i Vtsr'^, where V is the velocity there from right to left. 

1. Fixed Spherical Obstacle. 
We may assume 

t// = A + B cos ^ + (Cr^ + ^j I3 (cos 6) + ^, I3 (cos 6) +, &c., 

the centre being the origin. 

Then, if a is the radius, B = E 2= . . * = 0, and 

Ca^ + °- = 0. 
• a 



m'"^ 



But 

f %(cos 6) ^ — ^^ 6 



490 



MR. R. A. SAMPSON ON STOKES'S CURRENT FUNCTION. 



Hence 



xli = Y('i 



o 



1 

2 



a^ 



r 







2 (cos 0) 



a' 



sin^ ^ 



(58), 



2, Spheroidal Obstacle. 
Assume 

^^:-A + [BTo(p) + CH„(p)]Io(<z) 

+ [DTj(p) + EH,(^.)]iar/) 

+ [KIg ( p) + LHg (p)] L^{q) + . . . . . . . 

Putting D, K, &c., zero, in accordance with p. 489, the first terms become 



(59). 



smce 



B —- Cp + Hq, 
Eo{p)=p. Hj(j?) = l. 



Of these the term Cp must be excluded ; for referring to Chapter TV., p. 478, we see 
that the velocities normal to the surfaces p, q, are 



Q 



1 


/l- 


- (M'' 

I 


h^ ''^ 


¥ f - 


- q^ dq 


1 

—, M 


/I- 


- ]r d^ 



h 



xs 



r 



f dp 



e « 



(51). 



Hence the presence of this term would make Q infinite at the axis^ where 



Vi 



ih.^{l -f. 1 - (f) 



These expressions, in fact, tell us that where we may have 1 — g^^ = 0, the terms 
of t// that involve p, must contain the factor 1 — q^, and vice versa. The loci to 
which these values correspond are given on p. 478. They are, namely, for 

Ovary spheroid and hyperholoid of two sheets ;•— 

j9 =^ 1 on axis ; oo at infinity. 

g' =2 1 on axis ; at plane of symmetry* 

Planetary spheroid and hyper^boloid of one sheet ;— 

p =:z within ring of foci ; i oo at infijiity. 

q =:= 1 on axis ; on plane of symmetry beyond focal circle. 
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From these we see, too, that p is the member that becomes infinite, and hence, 
whether it occurs in the term In{p) Xi (?) ov not, or whether i/^ is a solution of Dt/; == 
or not, no higher positive power of p than the second can occur in i/f. 

Thus, in the region outside a spheroidal obstacle, we may have 



V* - A + Eri + [FI3 {p) + GH3 (p)] I, {q) + LH3 ( jp) I3 {q) + . . 



givmg 



17' 



L 



9 



0, 



^\,= 2Yhn,{q).\{p,) 



%(V) 



(P) 



I2 (I'd) H2 (j3o)_ 



. . (60), 



where P(^ is the parametei- of tlie spheroidal obstacle. 

When Pq is very small, that is to say, when the obstacle approximates to a disc, 



^-2VA^l2(^i) h{p) + 






a form which, shows that the velocity becomes infinite at the edge of a disc. 

A solution will be given hereafter, namely, on the supposition that the liquid is 
viscous, and the motion slow, which makes the velocity zero at the edge. See p. 504. 



3. Within a Ilyperboloid, 

Let us consider the motion within a hyperboloid of one sheet, q^, let us say. The 
velocity cannot, of course, be finite at infinity, and V will here denote its value at the 
centre of the hyperboloid. 

We may assume 

xlj = A+Eq + G.H, (p) I3 (q) + LH3 (p) I3 (q) . . . , 



and we find that all the coefficients must vanish except A and E ; giving 



xjj ^ Yh'q 



(61). 



It is easy to see that if we attempt to proceed to the case of flux through a circular 
hole in a wall, we find an infinite velocity at the edge. See p. 509. 



4. The Obstacle a Tore, 



The solution in this case, and in terms of the functions which I have been, dis- 
cussing in this paper, has been given by Mr. HiOKS ('Phil. Trans.,' 1884), who also 



3 K 2 
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considers the cyclic motion which can take place. No useful purpose would be served 
by quoting here Mr. HiCKs's results. I shall, therefore, confine myself to proving 
rigorously a point which he assumes, namely, that in the expression of the current 
function due to the motion of a tore through a liquid at rest at infinity, in the form 
(57), the functions H do not enter. 

This is deserving of some notice, since it is no longer true when the liquid at 
infinity is in motion. 

Assuming (p. 482) 

1 ?^ = 00 

'^ = Y^r — ^. tv ^ {^>^n+i {x) + B,H,+j {x)} COS {n^ + a„) . . (57), 

the condition that the velocity should be continuous requires dxfj/d^ =z at the axis, 
i.e., where x = I, 

Now I;^ + i(l) = 0. This condition, therefore, gives us no information as to how 
the functions I enter the expression. But it shows that if 

X = SB,,H,,+x (x) cos {n^ + a,,), 
then 

2(1 -cos£)^-Xsi^^===0. 
when 

for all values of f. 

In the first place, omit a^^. 
We find 

2(1 - cosf).0 - sinfBoH. 

— 2(1 — cos ^) BjHa sin ^ — sin ^B^H| cos f 

4- 
1 • • • 

— ■ 2(1 — cos ^) BJl^ + rU sin n^ — sin ^B^H,,+^ cos n^ 
+ &c., ad, inf. = 0, 

the alignment 1 being understood in the functions H. 
Hence equating to zero the coefficient of sin n^, 

{2n - 3) B,^iH,„. ~ inBMn^^ + {2n + 3)B,^iH,^. = 0, 

with one exception, viz., 

- 2BoH^ - 4B1H3 + 5B3H. = 0. 

But 

1 
H^^ + i(l) = — ^ 3 I ........ (45b), 
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whence 



2^ + 1 2n-"l 3 

Whence B^^ = 2Bq + i^^^? where Bq, yS are arbitrary. 

Similarly, if we supposed x to consist of terms B^^jH,^+j(^)sin n^^, we should find 
B'o = 0, B\ = 0, and, therefore, the whole series B'^^, = 0. 

Hence, if the functions H occur in xjj, they must be in the form of the series 

1 f 

a [I Hj {x) + H« (x) cos i + H| (x) cos 2| + . . .] 



(.2? — COS ^)* 

+ ft [H| (x) cos I + 2^H|. (^) cos 2^ + • • •] f 



that is, by p. 484, 



TT 

4 



The constant term, a — 7r/3/4, is irrelevant; the term 77/3/4^^ implies a finite 
velocity at infinity, viz., 7rft/2, from left to right. 



Chapter VI. — Rotational Motion. 

Among the problems on rotational motion in liquids, those are probably of most 
interest in which the circumstances most closely resemble nature, those, namely, in 
which the spin is due to the internal friction of the liquid. Certain cases of rotational 
motion in perfect fluids have been discussed by Mr. Hicks ('Phil. Trans.,' 1885), in 
connection with the theory of the motion of vortex rings, in terms of the functions 
I have before alluded to (p. 473). I confine myself, in the present chapter, to the 
consideration of certain cases of motion in viscous liquid. 

When the motion in a viscous liquid is slow, the equations may be reduced to the 
forms, (Basset, ' Hydrodynamics,' M73) 

dP oo dfj o d^ „y 

dt ^ ' clt ' dt 

where f, tj, i are the components of spin parallel to fixed axes, and v = fi/p, the 
kinematic coefficient of viscosity. If the motion is symmetrical about the axis of z, 
we may put C == 0, i = — sin <{>, tj = cd cos (f>, where co = — l/2t*r Dxfj is the 
resultant spin, and ^ is an azimuthal angle measured about the axis of z. But we 
have seen, in the introduction (p. 452), that 

'"'"^DV^: V^'^V (4); 



m w 
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write 2otco for V, and we get 



^ V at ^ 



Stokes's well-known equation, ('Cambridge Pliil Trans,/ vol. 9, Pt. 2.) 
Hence when the motion is steady, t|/ is a solution of 



Or we may assume 
where 



J}^ = . 



t/;=r V] 



» o ir Q tt -a 



• 00 



• 00© 



. (3 c/). 



(3e) 



so that t// consists of a function of the form of V, together with a particular integral 
of (3e). 

When the boundaiy is a symmetrical obstacle, the conditions which must be 
fulfilled there are : — 

(1) The velocity normal to the obstacle shall be zero. 

(2) The tangential velocity shall be proportional to the stress in the same 

direction, per unit area, across the surface of the obstacle. 




Now if u^ v^ IV are the velocities at any point parallel to fixed rectangular axes, 
whose origin is at the point, and if P, Q, R, S, T, U denote the stresses per unit area 
in tlie fluid, at the point, in the usual manner, as illustrated by the accompanying 
figure, then, it is known (see Basset, ' Hydrodynamics,' vol. 2, p. 242) 



du 



P = — . p — . ^[mO + 2ft — zzz — jp »- 1^^ -|- 2/xe 



d(X) 



Q = — . p — ^^0 4- 2/jl 



dv 
dy 



p-il^0+ 2/x/ 
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where 



clvj 



IE = - i? - f/.^ + 2^-- = - |) ■-- l/x^ + 2/^,9' 



clz 



^ /dvj dv 

\dy dz 

^ (du dw 

T = ^ (, rf, + 






^="(1+ 






2/xa 



= 2/xZ> 



2jac. 



„ c^?i- <:^y (i-z^ 

C7 ZZ= — -j~ -— -(- 



dx dy 



dyZ 



(62). 



I will now write these in terms of a, ^, y, the parameters of three orthogonal 
surfaces, whose intersection may be a convenient determination of the point. 
Let 

dx^ + dy^ + dz^ = A^da^ + B^d/3' + G'dy^ 




Then, if 0^, ^3, 0^ are the angles through which the normals at a, 0, y to the three 
surfaces must be rotated about themselves to make them parallel to the normals at 
a + doL^ /3 + d^i y + ^y? it is not hard to see 



B A 

do = .7- d0 — 77-77, da 
^ Adoi ^ Bd/3 



(63), 



with corresponding expressions for 61, 6^. For, in passing from P to h, along tlie 
normal to the surface yS, the normal to the surface a is rotated through the small 



d^ 



angle ebd^ or -rj dfi^ and, in passing from b to (i, it is rotated through the angle ead 



clA 



or :^Tr^ c?a, in the negative sense. 
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But the cliang3 in u, the velocity parallel to the normal at a, ^, y to the surface a, 
as we pass from a, /3, y to a + da, /3 + ^/3, y + dy, is 



du J dii 7^ du J 

- da + tt: <^p + v~ <^y 

doL did d^ ' 



vQo^- %oQ. 



2 



= d 



a 



du V dA w dA 



+ d/3 



die 



Hence we get 



e 



1^ 
A 



2a 



dtc V dA w dA 



v__d^' 
A dot 



+ dy 



'die IV dC 
dy A d 



1 dv 1 dia w dG v dB 



• • • 



Oi 



■ (64), 



with four others, which may be written down from symmetry. 

These enable us to express the second boundary condition (p. 494) in terms of any 
coordinates which we may wish to employ. 

Let us consider the solution of (3e) in terms of the coordinates r, 6. 

We may assume, p. 489, 

n — 00 



V 



n ~ 



and 



/73 1 _ ..3 ^2 

D = ~ + ~ ^ 



dr^ 



dfj} 



• « 



* • • • 



where /x = cos 9, 

Hence we require a particular integral of the equation. 



d^^^n ^^ (^ — 1) 



dr'^ 






i/;, == A,7- + B,r~^^-^i. 



Clearly 



^n = 



A,,r^ + s - Bnr-'' + ^ 



2 {2n - 1) 
Is a solution. 

Hence, a particular integral of (3^) is, changing the arbitrary constants. 

tl = GO 

^ji= t (A„r»+^ + B„r-^«+3) I„ (cos 6), 

n = 

and, omitting terms, which, by the expressions (p. 453)^ 



(6), 



H = 



1 # 

r^ Bill d do 

1 d-^' 

r Bin 6 dr 
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may be seen to give infinite velocities, we may write as a value of \jj, capable of 
expressing any slow, steady motion in a region of viscous liquid from which the 
origin is excluded, 

xjj = Di cos e + (C^r^ + D^r-i + B^r) I^ (cos d) 

9t = CO 

+ t (D„r-«-^i + B«r-«+8)L(cos0 .... (65). 



n = Z 



An expression equivalent to this is given by Mr. Butcher ('London Math. Soc. 
Proc.,' vol. 8). 

In terms of the coordinates r, 6, I shall consider two cases of motion, viz., (l) 
when the obstacle is a sphere, (2) when it departs little from sphericity. When 
there is no slipping, the boundary conditions are B; = 0, © = ; when there is 
sHpping the conditions that hold at the boundary are not easily expressed in the 
second case ; in its discussion I shall, therefore, confine myself to the supposition 
that the velocity at the boundary is zero. 

1. The Obstacle a Sphere, 
We find from (64) 

^ dn cm ® 

2c = — + ~ 

rdO dr r 

and, consequently, if a be the radius, we must have, when r = a, 



or, written in terms of i/^. 



R= 0, 

(dR .d% 0\ ^^ 



^=0 

de 



^ ^ ^ dfj? dr^ dr fi dr 

where y8 is a constant, which is infinite if there is no slipping, and /x = cos 0, 
But 

"^ ^^ ^" /^') dp? -^^ (^^ = n (^ - 1) I. (/x). 

Hence the above conditions require i)^ = B^ = 0, except when n = 2, in which case 
we have 

C^a^ + D^a-^^ + Bcfi—O, 

2C,a^-.D,a-"i(l+^)+B,a=0. 

MBCCCXCI. — A. 3 S 
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These give 



t|r = C 



a 



3(1 + ^f^\ ^ 



^s 



a 



% 






Cf 



' (> + i 



g (cos c?) 



66). 



It is clear that 0^ = V, the velocity from right to left at infinity. 

Passing from this well-known result, which is here given for reference and illustra- 
tion, let lis consider the motion about an approximately spherical, symmetrical 
obstacle, at whose surface no slipping takes place. 

j^Qj.^ ^ ^^ ^^l ^j'^Qjj \yQ ^j^0 equ:ation of the surface, and let f{0) be developed by 
Chapter IL into the form SaJL^^ (cos 6), 

Let us consider the surface 

r ^:^ a [1 4" ^n \ (cos #)], 

where we shall suppose a^^ so small that its square may be neglected. 



We must have 



d'\^ 



0, 



d^ff^ 



in (65) when r = a [1 + a,, I,, (cos 0)] ; and since the 



dr ' dd 

motion will not be far different from that in the case when the surface is a perfect 
sphere, all the cofficients which occur are of the order of a^^, excepting C^, Dg, B3, and, 
therefore, except where these coefficients enter, we may disregard the departure from 
a spherical form. We find, omitting the argument cos 6, and writing D'^^, B'^,, for 
D^a""'*"^^, B^a*"^"^^ respectively, 



a 



d'\}r 
dr 



0. 



j2C^a^ — D^a"-! + B^a) + a„, 1,, {2C^a^ + 20^^)] I^ 

+ S [ - (r - 1) DV - (r ^ 3) B'J I 



d^jr 



d (cos 0) 



C(C#^ 4" I)^a""^ + Bga) + ^fiXi {2C^a^ — D^a""^ + B^a)] cos 

+ 1 [d; + b;] p,. 



0. 



TT 

Jtlence 






203^^ 



Making use of these results we find 



whence 



for all values of r, and 



s[d; + b;3p„„, = o 



u f^ 4" t3 f •— ~ 



Ot^^ J.^ ZiXyiyCv Xa "~r* ^ ZiS3 f jLf «<—— u» 
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But 



y T ^n -r t ^n T f ^^ T 

A# % — — l^^+2 -r 5 hi "T" "5" J-#-s 



where S,,, €,,, 4 are given by (17), p. 460. 
Therefore 



tt 



i^ n-\-2 — J^2^ ^n n ? 






■/-i 



jDow 06;^ 



-» 



2 



Jj ^/ „_ o — — JDoCv 



a, 



b« 



3^ ^i^ O > 



all the others vanishing. 
Hence 



Ve^s 



Ig (cos ^) 



3r 
a 



2 

"~"^5 ^'^ Q ^« — 3 \^COS i/j 









a' 



a; 



n - S" 



' ■"" & /vjW O 



3Va 



3 



2 '^^ 



:;* 



2 



I^ (cos ^) 



a 



f^r 



W — 3 ^fi — 1 



3Va2 g, 



5t 



^?? 9 -^?H-2 (cos C7J 



tv 



w — 1 



,*j — 3 /i^Ji ~ 1 



a 



fi + r 



,^.n — 1 



w + 1 



r 



• • 



. . . (67) 



From this formula the solution for the case r = a {1 + totji I^ (cos 0)} can be 
built up. 

2. Spheroids. 

I will now consider the solution of (3e), p. 494, in terms of co-ordinates, p, q, 
defined on p. 477. 

By p. 479^ we may write 



% =: C» 



V=: 2 



CliXtvA 



D 



n = 



- 1 



A J,, (p) + BMn (p) 



Ifi. [q), 



Hence 



h^ (p^ — q^) 



■dp #" 



• • • • e 



(6a). 



(1 _ p2) ^V __ (1 _ ^2) f^=„ 

^ ^ ^ dp^ ^ -^ ^ d^ 

But, when n is not less than 4, 



'rt = 00 



h? t ip^-q^] 



n = 



A J^ (p) + BJIn (p) 



Jn (q) . (68). 



P 



s 



by pp. 460, 469. 



ln(p) 

H« (p) 






•t»-»-3 (P) 



O S -«y 
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Henc 



e 






I^ (P) 1.^ + 3 (2)] + ^nCn [I/i-S (P) In (q) — In (p) ln-2 (?)]• 



The corresponding terms in a particular integral of (68) are^ by p. 480, 



/V J\.f)G; 



n^n 



2 (2n + 1) 



[1^ + 2 {p) M^) + 1^ (p) I- + 2 (?)] 



2{2n- 3) 



[I« (i>) I«-2 (2) + l«-.3 (p) I. (gV 



Hence the terms arising in a particular integral of (68) from the terms on the right 
for which n is not less than 4, may be written 



where 



5t 



A^ 



» + a 



2 (2% + 1) 



\_^ffin ^n + 2 4^ + sj 



A2 



« + 2 



2 (2n + 1) 



\J^n^n J^n^2 4^4-2j' 



The terms on the right of (68), when n is less than 4, are 

ip^ - q^) . { Ao + B,p + [A,j> + B J g + [A,I, (IP) + B,H, (p)] I, (^) 

+ [A3l3(p) + B3H3(p)]l3(2)} 
= 2Aol2(p)+2Bol3(i>)-2[Ao+Bop]I,(g)+2[Ail3(p)+BA(i>)j2-2[Aij)+BJIs(^) 

+ A3S3 [I. (P) I. (?) - I. (P) I. (?)] + 63(33 [H,(i>)l3(2) - H, {p) I, (5)] - f I, {q)} 
+ A3S3 [I, {p) I3 (g) - 13 ip) I, (g)] + B3 {S3 [H, {p) I3 (^) - H3 (p) I, {qj] - -,\ . 13(g)}. 

Thus, the series of terms we have already found to occur in the particular integral, 
of (68), continues regularly down to ^ = 2, while there are also found the irregular 
terms 

- ^ol2 (P) -^ "^ I3 (P) 



II I I"! I ' 



HWili^H 



+ 



A, 



B1I3 ip) - IT h (P) 



K + [^o + j)p i^iq) 



q 



B, . BA . A 



v + 



45 



; + jp 



I3 (q)- 
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Thus xjj consists of these terms, together with the series 



n = CO 



t {9t„ + 3 [I«+a (p) In iq) + In (p) In+^iq)] + S^+s [H„ + 3 (i?) In {q) + H« (p)I«+a(?)]}, 
and the series 

H = 00 



n = 



All the constants above are independent. But referring to p. 490, we see that the 
conditions for finite and continuous velocity everywhere require that Aq, Bq, A^^, B^ ; 
3I45 5(5, &c.; Cq, Dq^ Og, C45 &c., should all vanish. Hence all legitimate expressions of 
\fj are included in 

t// == Co + Diq 

+ I2 (q) [B^P + 0,1, (p) + D,H, (p) + S,H, (p)J 
+ I3 (q) [B3 + D3H3 (p) + ®,H3 (i.)] 

+ I4 (2) [® A (i^) + I^^H, (i)) + «eHe (i>)] 

+ &c., regularly (69), 

writing B,, B3 in place of -^, t| respectively. 

The first case to which I shall apply this formula is that in which an obstacle in 
the shape of a spheroid is opposed to the flow of liquid, and, to begin with, I shall 
suppose there is no slipping at the surface ; so that the boundary conditions are 

^==0 ^=0 

dp dq ' 

when p =Po, where p^ is the parameter of the fixed surface. 
These require, omitting the argument, p^, 

B,p, + G,l, + D,B, + 93,H, = = B^ + C,V, + B.Q, + ^.Q, 
Bg + D3H3 + ^gHs = = D3Q, + aSsQ, 



Hence, three of the above constants are arbitrary, namely, we may give what 
values we please to C,, ©4, Bg, and all the others will then be determinate. If we 
give a finite value to SS^ or B3 we get an infinite series which in each case is diver- 
gent (see p. 508). Let us then first consider the solution given by aS^ = 0, B3 =: 0. 
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We have 



where 



Hence, since 



we get 



where 



^ = h (?) [CJa (p) 4- D3H3 (p) + B^p] 

C#o+ D3 Qi (|)o) + B3 = 
Csia (Po) + D2H3 (po) + Ba^o = 0. 






-Po 






3 



T> . 

X>3 — 



2L (po) ^'^' 



c< 



2L (po) 






Now the term which gives a finite velocity at infinity is Cgl^ (p) I3 (q) = •-■ Q^m^jih^, 



Hence, if this velocity be V from right to left, Cg = 'iWY, and 



V. = 2 w I, (2) 



ijp) 



i^ + ip^' + 1) H, jp) 
2L (i'o) 



• » • • » 



(70). 



From this result we find, by simple reductions, that the velocities U, W (p, 478), 
perpendicular to the axis and parallel to it, are 



U 



V g (2^^ - i^o^) \/(l - <f) 



w 



Y 



2L (j?o) 



' + ^Moa'^--+ U ^^f^> 2L (po) 



2 



i^ 



1 



i^ — g; 



(71), 



From (70), by the equation t/; == const., we get the forms of the stream surfaces, and 
from this point of view, an integral which includes (70) has been given by Mr. Herman, 
' Quarterly Journal of Mathematics,' ] 889, No. 92. Taking Oberbeck's values for the 
velocities due to the steady motion of an ellipsoid through a viscous liquid (' Crelle,' 
vol. 81), he finds as one surface on which the stream-lines relative to the ellipsoid lie, 



1 + 



d\ 



0^ + V {{a? + X) (6^ + X) {g^ -r X)}^ 



(A, + B, + C,) € 



xy = constant, 
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where 



1 00 



A, 



and (a^ + e)^ is the semi-major axis of the confocal to the moving ellipsoid, which 
passes through the point x, y, %. I have slightly changed his notation, making the 
axis of z that along which the ellipsoid is moving. 

But if the ellipsoid is one of revolution, so that a = ?>, then 

K^a = loff ^'^-] - - 

p — 1 p 



.00 



in the notation of p. 477, where 



Hence, writing Mr. Herman's integral 

[Ho + c^Co - H. - c^a -■ € (A. + B, + C.)] xy = const., 



and substituting, we find the expression in brackets equal to 



1 f 



4p 



p^ 



I |4L (po) - 2 (|>o^ + 1) [i log^ ^ ^ 

and 

xy = 73r^ sin ^ cos ^, 

= A^ (p^ — 1) (1 — q^) sin ^ cos <^, 

^ being the azimuthal angle about the axis. 
That is to say, the surfaces 

[2L (po) Is (P) - {Pd" + 1) Hg (p) - j9] I^ (g) sin (^ cos (^ = const., 

are composed of stream-lines ; and the motion being symmetrical about the axis, 
(j) = const, is the other family of surfaces composed of stream -lines. This gives the 
expression (70) equated to a constant, as representing stream surfaces. 

From the expression (70), we can obtain the expression (66), when yS is infinite, 
giving!/^ when the obstacle is a sphere. We must suppose h indefinitely small; then we 
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find p = - , g = cos 6, where 6 is the angle between the asymptotes of the hyperbola q, 
and is ultimately the vectorial angle of any point. But by (28), when h is small, 



and 



H, {p) = - 



h 



r 



L(^o) = H3(po)+ ilog- 



P + (i^o'+ 1)^3 (p) 






1 



r 



a-^ 



r\s = ;?/V sin^ ^ 



T — 



2. 
3 



/^ 



a 



a" 






a^Y sin^ 
4 



8^ ci 2r^' 



a r 






which is the expression of p. 498, in which /3 is supposed infinite. 

Three particular cases, besides the sphere, are worthy of notice. 

(1.) The obstacle a paraboloid^ of parameter 21^. Taking 2Z, 2m as the parameters 
of the confocal paraboloids through any point, the solution is found by making h very 
great in (70). 

We find, without difficulty. 



^ = 



V 



2 log Ahllfy 



m 



llog^j + l-l 







. • 



. . . (72). 



(2) A planetary spheroid of eccentricity l/y^2. Hence, by p; 478, 



Pq — ^, 



and (70) becomes 



rl, = a« V I, {q) 



^Ap) + :- 



t 

i 






e Q • 



where 2a is the greater axis. 



(73) 



(3) A circular disc. Here p^ = 0, 



Zi% 



h{p)-^ -{p + '^Av)} 



• 



(74), 



where a is the radius. 
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Turning to (51) we find the velocities given by this formula are 



U = - 



2V qp^ v/(l - q^) 



iir (j)^ — if)\/(p'^ ■— 1) 



w 



2V 



iw 



'^ ^ 2^ -—1 



pg- 



2r — g 



ITT 

~9. 



vanishing at the edge. 

As an example of the analysis to which the use of the expression leads, I will now 
find the value of xfj for a surface approximately spheroidal, say p =Pq+ a,,ln{q)^ 
where a,, is a small quantity whose square will be neglected. 

First suppose n an even number ; let us then assume 

'/' = I. (q) [B,P + CJ, (p) + D3H3 (p) + 33,H, (p)] 

-f- . . . 

+ l^ k) [^n'Hn^2 (P) + ^nBn (p) + SB^e^-sH^ + S (p)] 

at which point we shall find that we may stop ; in fact, there are n + 3 constants in 
the above expression, and n/2 + 1 terms, each of which yields two equations corre- 
sponding to d\fj/dp = 0, d^/dq = ; we are thus left with one constant arbitrary, which 
is determined when we know the velocity at infinity. It is clear that all the constants 
except Bg, C^, D^ are of the order of a^, and hence, at the surface, we may put p z= p^ 
in the terms where they arise as coefiicients. 

The condition -j- = when p =i p^ -^ a^ln (q) gives 

1 [B2P0 + C3T3 + D3H3 + 93^H J 

+ a„ql,, (q) [B, + D,P, + D,Q J 
+ ^B {9) [^4^3 + D4H4 + ©eHe] 

I • • « 



dp 



gives 



MDCCOXOI. — A. 



I2 (2) [B, + C,Pi + D A + a3,Q3] 

J ... 

+ L + 3(?) [®/^ + 3Q^^-l+ D,,+3Q,, + J = 0, 

3 T 
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Now, 



i.(^)i.(?) = %i».s(<z) + V^i42) + |i»-ste) • . . . (17), 



9 — «-t-a \a:/ i o 

by p. 460. 

Hence the second series of equations requires 



B, + C^Pj + D^Q, + a3,Q3 =■- 0, 



« • 



a • e 



®^* + 2^^^- 1 H" MiQ« + 1 + ^ Oft 



0, 
0, 
0, 



where 



a. 



2 



C. + D 



fZQi(^oy 



2 



C^jPi 



J 



Making use of the first of these^ we see that the former series requires 

^^0 + C,I, + D,H, + g3,H, = 0, 
©,H, + D,H, + SeHe - 0. 



« 9 



'^n + 2-*^« I ^^i + 2-*^-^' + 2 = U , 

Taking any pair of equations 

®rH^-2 + D^H,. + 3SrH.2Hf+2 ~ 0, 

where r Hes between 4 and ^2 — 4 inclusive, we find 



where 



D^.V^. = 33^.U^. = 93r+3U^-f-2 = ^85 a constant^ say, 

U ^- = Xl ^\olr _ 3 "~" ±l.r _ 2 v^r _ x j 



In place of these, the last three give : — 
Last pair, 



'L'w + 2^?« + 2 — ^^« 



H 



J-*^« + 4 ' 



M + 3 



^^ + 2 



U 



» + 4 
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Last but one, 






U^y ^ — — '<Jft^n ^ \^n """ ■*-/ 



-« 4. 3 * « 



u 



w + 3 



Last but two. 



'n^^ n 



XX ■IV Y 4, 



I)u^^Vn-2 = /3+ OtL 



-?i. » « — 3 



IJ 



We have here two values for 33,^ + ^U^^ .,. ^ ; reconciling them, we find 



or, 



= ^ + al,JIn^2 + Gt{en — l) H« + aS^H^^ + g, 



With regard to the first pair of terms, we suppose Cg arbitrary, and Bg, D^ 
separated into two parts, of which the former are the same as when the spheroid 
is perfect, while the latter are of the order a, and fall in regularly with the above. 

Hence we get 

^ = I,(g)[Bp+CI,(p) + DH,(p)] 

P 



+ 14 (2) 



jj,^ V, ■+■ u, . 
H3 (p) n^jp) H(, Oj)' 



u 



4, 



V. 



u, 



"y" * 



TJ * V 



^"^ 1?« — 2 \2/ ^^^ IT -ti^i — g^^/ b?i TT Xlfi \pj 



+ 1. (9) 



u, 

tRi-i (Po) 



n 



u 



H«_2(p) 



n 



+ ( e, _ 1 ?^iM _ ^-"-^^^o)^H4p) ^ ^^^M H,^^ (p)' 



+ 1^ + 3 (?) 



u 



« + 2 



V 



u. 



+ 3 






In this expression, B, C, D have the same values as the corresponding coefficients 
in (70), and 



3 T 2 
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If n is odd we get a series exactly similar, except that the first terms are 



« W - 1) H« (Po) 1 13 (?) 



U3 "^ V3 "^ Us J 



+ &c. 



There is one exceptional case, viz., tc = 1 



We have 



where 



V.= r,(2)[Bp + CI,(p) + DH,(^3)] 



+ 



u 



Q3 to) 



[Hg (i5) - H3 (po)] I3 (^) , 



a = — ex. 



C + D 



c* 4. T) '•^ ^i (yp ) 



:/^o 



I shall now consider those two other functions^ which, as we saw on p. 501, are 
solutions of J)^ = 0, and make -^ = -^ = 0, when p = Po- 

One of them consists of a series of terms of odd order only, and the other of terms 
of even order only. 

Consider the latter ; put 02=63 = 0. Then we have the unlimited series of 
equations 

Bapo + D,H, + a3,H, = =B, + D,Qi + ®A 
95,H, + D4H, + aSgHfi = = SB4Q1 + D4Q3 + SScQs ; 



C • • • e a 



whence, as on p. 506 



■^/i ^ n — '-^^n ^ n — '-^n 4. g U ;^ ..j. g 



B,U,. 



And we find 



'/' = 



B3U, 



+ 



_^ ■ Hg {p) H^ (ff) " 

lU,^ V, "1" U. 



4 



H 



TT •" V TT 



I2 (?) + 



'^ + 3 



I;, (g) + &C.J ad inf. \ , . (76). 



Similarly, the series consisting of terms of odd order is 



^ 



B3U3 



1 jL. IsM . liM' 

lUs"*" V3 "f" U, . 



^.3(2) + '^'^■' '^'^ ^-^i*^ r • • (''6a) 
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Now, the maximum value of the term 

H,, (p) Mn (p) 



U« H,, (po) Q,, ^ 3 (pq) - H,, ^ 2 (pq) Q«_i (pq) 

occurs when p zzip^; and this maximum becomes very great when n is very great ; 
for by p. 475, when n is great H,^ (x) = — |- a/ --^ I^'*""^ (1 — i^f, and (Heine, § 40) 

where io = Po -— %/{Po^ -- i)^ mod. i^^ being less than 1. Hence this expression 
becomes infinite with n ; and no infinite series in which such a term is present can 
represent a motion within the limits of our discussion. 

We could have foreseen that these series would be inapplicable, for we have already 
found one to which there is no objection, and it has been proved by Helmholtz 
(^ Wiss. Abh.,' vol. 1, XII.) that there is only one slow motion which satisfies the 
equations for a viscous liquid, and obeys definite boundary conditions. 

3. Jtiyperbotoids. 

Let us now consider the motion which may take place within a hyper boloid of one 
sheet ; say, of parameter g^. 

The conditions dxfj/dp = d^jdq = are evidently satisfied by taking 

t|r = Dj^g 4" B3I3 i<l) 
where 

Di + B3P, (go) = , 
or 

and 

B3 = - 2VA73 (1 - <7o^) {n) 

That is to say, the stream surfaces are hyperboloids of the system confocal with the 
boundary. 

The velocities which {77) gives are 



"^ - l-g,' (p^ - cf) i V 1 

TXT _Y \T 9o) 

whei*e V is the velocity at the centre from right to left. 






J 



. . . . . (77a) 
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In the particular case of a flat wall pierced by a circular hole, we have q^ = 0, and 



U= - V-7 



n^ \/(i - i^) 



ii^f- (f)^{\ — f) 



W = V 



both of which vanish at the edge. 



II 



f -t 



(77b) 



4. Slipping Motions at the Boundary. 

Approaching now as near as may be possible to a discussion of the solutions when 
the liquid slips at the surface of the obstacle, we find, without difficult or very long 
calculation, that if on p. 495 we put a, jS, y, = p, q, <j6, and consequently u, v, = P, Q 
of p. 478, while w = 0, we find, I say, a = & — 0, 



2c = 



{f - i') V{f -1.1 -fY 



(i-3^)^-(p^ 






2g (1 - (f) df 2p(p^-~'Ji) d^ 
p^ — gp dq p^ ■— ^ dp 



(78) 



Now if T is the tangential stress per unit area at any point of the obstacle, whether 
spheroid or hyper boloid, T = 2/xc, and a boundary condition being that this tangential 
stress shall be proportional to the velocity of the liquid, at the surface of the obstacle, 
we have 

for spheroid [p^) T = ^Q, P =: 0, when p == p^ 

for hyperboloid {q^) T = ^P, Q = 0, when q = q^ 

where ^ is a constant, becoming very large when the slipping is very slight. 
Now 

rP=C,+ J),q + I, (q) [B,p + C,I, (p) + D,H, (p) + ^,H, (p)] 

+ I,{q)[B, + I),-H,{p) + ^,-H,{p)] 

+ I, (q) [aS.H, (p) + D,H, (p) + 33eHe {p)]+ . . . (69) 



and remembering 



dx 



• ^u \*^) — ^n^\ y^J 



fP Xfi (x) qi {^qi — 1) hi \^) a 

we see that if P, and, therefore, d\jj/dq = 0, for definite values of p, we also have 

(Pxfj/dq^ = 0. 

Hence the above conditions become, for spheroidal obstacle, p = p^, 
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dyjr 

dq 


= 0. 


C^^i|r 


M 


dy\t 


/3 


clp^ 


Po' - 


q^ dp 








d-xjr 

p^ — 1 diJ 

when 

2^ = Po5 (79) 

/or hyperboloid, q = qo , 

From these equations we can detect at once a solution in a particular case, viz., 
the flux through a circular hole in a plane wall, when the result we have already 
obtained on p. 509, viz., \// = — Va^S . q^ {a being the radius of the hole, and V the 
velocity at its centre) still applies. In fact, the tangential stress at any point of the 
wall is zero. 

But it is obvious that it is out of the question to apply, in general, the direct method 
of assuming (69) and substituting, in order to find functions which satisfy (79) and 
(79a), and it is also clear that the solution for the spheroid is not, as in the case of the 
sphere, of like simplicity with the case when the slipping is zero. In one case I 
have followed out the process of solution, viz., when the obstacle is a circular disc, and 
the slipping slight. We then have Pq == 0, and 

dsfr fi d^yjr 
^ dp /3 dcf 

and assuming that the solution does not differ greatly from the case when there is no 
slipping, I at first neglected, on the right, all terms except those entering in (70), p. 502. 
This led to a set of equations for finding the coefficients 

«.U. ^ ©.U. = ~ 5 . t . ^^""'^ 



'5^5 -^3^3 — *^'3 



33,U, - 93,U5 = - -V . (1)^ 



/3 



/3 



m TT m TT — (^ - •^) (2^* - 2) /2«, + 4. . . 2\3 2/xa2V 

Now, when n is very great, the limit of 



Hence, when n is great. 



2 . 4 . . . 2>i, . ^ I .IT 



33 u - ^ U ~ - ?-'^-- 

P 
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That is, the coefficients 33 become infinite when n is infinite, and the series which 
the expression (69) would require us to use for ^ is divergent, showing that the 
assumption that \jj may be written in a series proceeding like (69) is, in this case, not 
justified. 

It will be seen from the last section, and from p. 508, that an infinite series of the 
form, satisfying common boundary conditions, frequently becomes divergent. This 
would not always seem to show that there is no corresponding possible motion, but 
rather that the mode of expression chosen is not permissible. 

Now it is easy to see that there is no solution of D^^ = 0, consisting of a finite 
number of terms of (69) which obeys (79) or (79a), the solution \}j = q^ excepted. 
Hence it seems probable that a series of the form (69) is not appropriate for expressing 
the motion of a viscous liquid, slipping at the surface of a spheroid or hy perboloid. "^'^ 

Note. (Added, June, 1891.) 

At the suggestion of the referees, I add a more detailed consideration of the motion 
past a disc, and through a hole in a wall. 

A circular disc, of radius tf, is moved perpendicular to its plane with velocity V; 
the stream-lines relative to it are given by 



'/' = - Salvia (q) 
or, writing p = ip, 



I.(P) + ?{P + H,(|>)}' 



TT 






w 



The stresses at any point of the liquid are 

P=— P + 2/X6; Q=— p + 2/x/; R=-~p + 2^^; S=0 = T; U = 2/>tc, 

and p being the mean of the pressures normal to any three mutually orthogonal 
surfaces, I shall hereafter allude to it as the '' mean pressure '' at the point. 
Let us first find p. 

If Q = — . ^ — V, we have t 

P 

where ^ = \)^. 

Transforming to the elliptic coordinates p, g. 



do = — 7 
h 



—■ ([ d'p ^ 1 — ^r dq ^ \ 



* See p. 515. 

t Basset, vol. 2, p. 262. 
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which is an exact differential by (Se) p. 494 ; and gives, in the case of a,ny spheroid, 

Q _ „J^^_ . £„„„ 

or in the case of the disc, of radius a, 

4:vY q 



Q- 



O 9 



and since there are no applied forces, the variable part of the mean pressure is 



4Vyw, q 

^ TTCC ^^ + q 

Again, by the formulse, p. 496, 



1 /p^ — 1 dtt V q \/ (1 ~ t) 



0\ 3 J 



h V j>9^ .-- q"^ d]i Ih (27^ — q") 
1 / 1 ~ q^ dv u p \/0^^ "" 1) 



f/ 



^9^ — 2^ dq h {p^ — q^y^ 
up vq 



■7\3 7 



0\ J 



-y/ Q:)^ — 1 . p^ — q^) v{ 1 — 2^ . p'^ — q^) 

a:=^ h ■= ; 
2 --. 1 . /l-Tl-t . ^ I 1 . / tjul . ;^ _ t:;^ \/(/ - 1) , iiq \/(l -- g^) . 

whence, in the case of the di^^c, 

_ 2y ^ fq (1 4- g^ - 2p') _ 2 V ^ ]r£ 0_+_£jiV) . 
^ "^ ftTT ' (/"- <?2)2 (i92 - 1) "^ «7r ' ("))2 ~-f 2^)2 (:p^ + 1) ' 



2_v ^ £f- /I _ £^ _ 2 y ^ ^^/- / 1 -g^ 



c 



These all vanish on the surface of tbe disc (p = 0). Hence the total pressure on 
the disc is got by integrating the mean pressure over both sides of its surface.- This 
gives 16V/xa, agreeing with the knowai value for any ellipsoid/''' 

I append tracings of the curves , const. = ^^jci^^ P/'^l^^ ^) /> 9^ ^> with numerical 
values calculated on the supposition that V/7ra ^ 1. The first shows the stream-lines ; 

^ Bahret, § 497. 
MDCCOXCr.-A. 3 XT 
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the second, the mean pressure at any point ; the third, fourth, and fifth, the departure 
from this mean in directions respectively normal to the ellipsoids, to the hyperboloids, 
and to planes through the axis ; the last shows the tangential stress for these 
particular directions of normal pressure. 

With the exception of g, which has a maximum at the point on the axis distant one 
radius from the disc, all increase if we approach the edge of the disc in particular 
directions, and we can make this increase very large ; but at the edge itself they 
appear to be indeterminate, not infinite, all the surfaces of equal stress passing 
through one and the same line; and there seems no ground to suppose that it is 
impossible to impress a sufficiently great external pressure upon the whole liquid to 
prevent the pressures becoming negative. But it suggests that the test of zero or 
finite velocity at a sharp edge is not sufiicient to ensure the motion being admissible. 

Next, consider motion through a circular hole in a wall. 

We have, for any hyperboloid, 

where V is the velocity at the centre. 
Hence 

and the mean pressure is 
and for a flat wall (g'o ^ 0) 



tan ^ :p + "-^ 



r + r^ 



AuN 



a 



^^ + T 



tan ^ ^ + -^ - 



where a is the radius of the hole. 

The difference between the pressures at infinity is 

4/xV7r 



ill (1 - gj) 
Again, if F is the total flux per unit time through the aperture, 

Hence, a being the radius of the aperture, so that a = ih ^/{l — ^q^) 



"P ^Jt . (^gp + 1) (1 - qj ^ jj 
12/^ (1 + 9.0)' 



or, in the case of the plane wall, 



a^ 
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By applying the formulae given at the begmning of this note 

^ _ TSfi p(f a -\- if - 2f) _ 2V i^ff (1 + cf -{- 2^^) ^ 



/ 



a (2f- ^ ff {f ^ 1) a (V' + q'Y {f + 1) ' 



^ '^ 7" ' (i>' - ?) (p' ~ 1) "~" a'(f + q^) {p^ + 1) ' 

Tracings of sections of these surfaces are given with numerical values on the sup- 
position that Y/a = 1, as also for j9/2/>c. With the exception of the mean pressure, 
the curves are in general character the same as for motion past a disc. 

The surfaces of mean pressure are the same, except for the numerical values 
attached to them, for all planetary spheroids and for all hyperboloidal boundaries. 

Since the reading of this paper, I have found a proofs that the expansion of any 
function of q in a linear series of I (qYs is convergent where q is less than unity. 
This throws more light on the problem of shpping motion past a disc, attacked in 
Chapter VI., for it shows that if the expression for any function i// is a divergent 
series, then there is no corresponding finite motion. 

The work there was conducted on the assumption that the motion, if existing, did 
not differ greatly when the slipping was small from the case when there was none ; 
and, therefore, we have not a rigorous proof that there is no such motion, though this 
seems to make it highly probable. As it is not easy to see how there should be a 
motion with slipping for one particular case of spheroid, viz., a sphere, and none for 
another, viz., a disc, unless the former is a wholly exceptional case, it is interesting 
to notice that the probability of it being exceptional is borne out by the analysis 
(p. 511), for if we approach the case of the sphere, by putting p = r/h where h is 
indefinitely small in the equation of condition (79), q ultimately disappears from the 
equation. 

It must be remembered that the search for steady motions in viscous liquids 
proceeds on an assumption, viz., that if the motion is started from rest, the limit of 
the rate of change of velocity is zero at all points when the time is infinite ; therefore, 
such motions must always be exceptional, and. again, the analysis bears us out in 
expecting this, for the same thing will be noticed in the two solutions which I give 
that Stokes pointed out in the case of the sphere — that the possibility of satisfying 

the boundary conditions -^^ = = ~- depends on the presence of an irregular term^ in 

the series (p. 501), namely, Y>^f\^{(q) for spheroids, and '^^^{(q) for hyperboloids. 

* Inserted in Chapter II. 
3 U 2 
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Fig. 1. 
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